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ABSTRACT 


Present work deals with the analysis of a capacitor - diode 
voltage rnultiplier. The full wave doubling rectifier has been 
chosen for this purpose. Analysis and simulation has been done for 
the rectifier fed by ideal and non ideal sine and square wave 
voltage sources. Analysis has also been done for the multiplier 
fed by ideal sine and square wave current sources. Later in the 
thesis, a modification has been suggested to the classical 
doubling rectifier circuit, in \diich the diodes have been replaced 
by switches. This has been done with a view to controlling the 
average output voltage. Also studied is the voltage doubling 
rectifier fed by a square wave voltage whose pulse width can be 
controlled.lt is found that the average output voltage is not very 
sensitive to the variations in the pulse width of the square wave 
input voltage. Another scheme in vhich the classical doubling 
rectifier fed from PWM current and voltage sources, has been 
stvdied. Simulation results have been presented. 
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INlKODUCriaN 


1.1 V(X.TftGE MOLTIPLIERS AND ITS APPLICATIONS - AN INTOODOCTION 


Capacitor-diode voltage multipliers (CDVM) have been iJised in 
several applications as a simple and reliable means to obtciin a 
high D.C voltage from an A.C source. Among the advantages offered 
by this solution is the fact that the a.c voltage may be 
svdbstantially lower than the desired d.c voltage and that the 
capacitors and the diode ratings are lower too, as the sharing of 
the voltage stresses is intrinsically ensured by the multiplier 
operation. 

When a transformer is used as an a.c source for the 
multiplier, substantial benefits in terms of voltage and 
insulation ratings of the transformer itself are obtcdned. 

The CDVM advantages become particularly a^reciable in those 
applications, such as space and comnrunicatioans, where small size 
and weight and high efficiency and reliability aure of prime 
importance [4]. In all these cases, the main disadvantage of CDVM, 
i.e., the delay between iiqnit and output and normegligible amount 
of capacitance needed, arereduced within acceptable limits by the 
increase of the operating frequency of today's converters. 

In resonant converters, vdiich are used more and more owing 
to the advantage of efficiency and the effects of parasitic 




parameters they offer, the auiopticn of a voltage multiplier with a 
rescaiant ccsnverter may turn out to be very ccmvenient. 

There are many applications that require a well stabilised 
and regulated HV or BHT sij®>ply. Table 1.1 lists some of the 
important applications. In addition to all those listed in the 
table, HV sv^jplies are used for image intensifier tubes, camera 
tubes. X-ray machines, insulatuion testers, and other specialised 
tubes and instnments. 

1.2 MBTBODB TO FRODDCE HIGH VOLTAGE O.C. 

a) To produce d.c output above IkV it is possible to use a mains 
transformer with a secondary wound with the necessary nixnber of 
turns, and this is the solution used in some oscilloscopes and 
X-ray systems where a simple fixed voltage unit at low cost is 
required. The drawbacks to this method are that the unit is 
rather bulky, the transformer and the smoothing capacitors 
required being very large, and the HV generated at the transformer 
secondary is highly dangerous since it has a relatively low 
inpeadence. Smoothing capacitors in such units, because of their 
relatively his^ value, nay also store a lethal charge for sometime 
after the mains is switched off. 

b) Voltage Multilpliers circuits, of lihich there are several 
types, are very useful for producing high values of d.c voltage at 
low currents. Typical voltage multipliers are shown in Fig. 1.1. 



Applications 


Typical 

Current 


Typical Typical 
Vol tage Regul at i on 


1 . 


Cathode Ray Tubes used 
in oscillosccqpes, radar 
displays and VIWs 


A few ndlli 2kV to 
ani>eres 20kV 


2 . 


Photo multiplier tubes 
used in scintillating 
counter, flying spot 
scanners and low level 
photometery 


0.5 mA to 
5.0 mA 


IkV to 
3kV 


3. 


Lamp supplies for 
photo ccqpiers 


SidA 


5kV to 


lOkV 


TABLE 1.1 


1 % 


0 . 1 % 


few 

percent 







The full vrave doubling rectifier circuit (Fig. 1.1a) is 
basically two half wave doubling rectifiers with the output summed 
by two capacitors. ccsnducts on the positive half cycle to 
charge C to E (the peak value of a.c secondary) and then on the 

1 m 

negative half cycle conducts charge to also to E^. As with 
other rectifiers of this type, the regulation is poor so that only 
a relatively light load of a few ndlli amperes can be si^plied. 
Another doubler Fig. 1.1(b) commonly used is also shown. 

In this circuit must be a larger capacitance value than 
C^. on negative half cycles, D conducts and C diarges to E,. Cto 
the next positive half cycle, the left plate of goes positive 
and conducts carrying to almost 2E^. laical ly is twice 
the value of to allow correct chaurge sharing between the two 
capacitors. A voltage tripler can be constructed by adding an 
extra half wave rectifier as shown in Fig. 1.1(c). Alternatively a 
tripler and quadrupler in Figl.l(d) can be made using a few 
sections of the well known Cockroft-Walton cascade rectifier 
Figl.l(e). Because the output iiipeadence of multipliers increases 
rapidly with the nxirber of multiplier stages used, in most cases 
the rectifier in HV circuits is limited to a doubling rectifier 
circuit. T!he output inpedance of the voltage multiplier increases 
approximately by the cube of the number of stages [5]. A quadrrpler 
has an output inpedance about eight times greater than a doubling 
rectifier. High valxies of the output inpedance limit the 
available output current and degrade the load regulation. 



HV power units are specified in the same way as other d.c 
power siflpplies. Parmeters such as load and line regulation, 
output inpedance, output ripple, ternperature coefficient and 
stability being the most inportant. 


In this thesis analysis and simulation has been done for a 
commonly used voltage multiplier circuit ,the full wave voltage 
doubling rectifier. 

In chapter 2, analysis and simulation of a voltage doubling 
rectifier fed by ideal and non ideal sine wave voltage sources has 
be^ done, Differential equations were written and solved 
analytically for the rectifier fed by ideal sine wave. 
Differential equations for the doubling rectifier fed by non ideal 
sine wave voltage source are of second order and hence solved 
numerically. The ntmerical solutions have been detained using a 
subroutine D02BBF from NRG library. 

In chapter 3, analysis and simulation of the voltage doubling 
rectifier fed by ideal and nonideal square wave voltage source has 
ben done. Analytical solutions of the differoutial equations have 
been obtained for the former cause. Numerical solutions have been 
obtadned for the latter case using , again, D02BRF subroutine from 
NAG library. 



In chapter 4, analysis and simulation of a doubling 
rectifier fed by ideal sine and square wave current sources hais 
been done. It has been found that the rectifier works as current 
halver. 


Chapter 5 is devoted to SPICE simulation. In this cahpter, 
SPICE simulation of voltage dovibling rectifier fed by a square 
wave whose pulse width can be controlled, has been done. This is 
done with a view to controlling the output voltage. Also studied 
is the voltage doubling rectifier fed by FWM voltage and current 
sources. A modification has been suggested to the claissical 
voltage doubling rectifier. 2Vs per this modification, the diodes 
have been replaced by switches. Hiis modified doubling rectifier is 
found to generate controlled output voltage vdien fed from a square 
wave voltage sorurce of a constant pulse width. 



CHAPTER 2 


VOLTAGE DOUBLING RECTIFIER FED B3f SINE WAVE VMiTACaE SOURCES 


In this chapter analysis and simulation of the voltage 
doubling rectifier fed by ideal and nonideal sinusoidal voltage 
sources has been done. Simulation results have been discussed. Some 
inportant coiclusions have been made. 


2.1 VOLTAGE DOUBLING RECTIFIER FED BY 

AN IDEAL SINEWAVE VCX.TAGE SOURCE 

The doubling rectifier (Fig. 2.1) vd.th ideal sinusoidal si;pply 
has three modes of operation. Hiese modes are characterised by 
the patterns of diode conduction. The mode transitions are shown 
in Fig. 2. 2. The modes are: 

Mode 1 :Di ON , Dz (XT, a< «t i /5 

Mode 2 :Di OFT*, Dz OFF, p < at i n + a and n + fi< at ^2n+a 

Mode 3 :Di OFF, D20N,n+a<<ot^n+p 

The description of each mode is eis follows: 

Mode 1: In this mode the diode Di conducts as svpply voltage is 
larger than capacitor voltage , Vci at ut = a . The equivalent 
circuit for this mode is shown in Fig. 2.3 . At wt = p , diode Di 

goes off as the soxirce current, i« becomes zero and circuit 

transits to Mode 2. Diode Dz always remains off throughout 
this mode as the voltage across its anode to cathode never becomes 
positive. The voltage across capacitor 1 is same as the supply 
voltage. 



Fig. 2. 1(a) 


Fig. 2. 1(b) 



Voltage Doubling Rectifier fed by an ideal sine wave 
voltage source. 



Circuit Waveform for the circuit of 

Fig. 2. 1(a). 


Mcxfle 2: In this mode both diodes Dt and D 2 are off and current 
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supplied by source is zero. Voltage from anode to cathode for both 
the diodes renains negative throughout this mode . At ut = n + a , 
the voltge across diode D 2 becomes positive and this diode turns 
on. The circuit transits to Mode 3. The equivalent circuit for Mode 
2 is given in Fig. 2. 4. Load is stjpplied by energy stored in 
capaciotors Ci and Cz. 

Mode 3: In this mode the diode Dz conducts as the negative of 
si:5>ply voltage is larger than capacitor voltage , Vcz at 
(i)t = n + a .The equivalent circuit for this mode is shown in 
Fig. 2. 5 . The voltage across capacitor 1 falls . The voltage 

across capacitor 2 is same as the negative of the si^jply 

voltage. At wt = n + p , diode Dz goes off for the source 

current, i* becomes ec[ual to zero and the circuit transits to Mode 
2 again .The diode Di remains off throughout this period from <>)t 
= n+ atojt + pas the voltage across its anode to cathode 
remains negative . The circuit remains in Mode 2 till {ot< 2»t + a 
and at wt = 2it + a the circuit transits to Mode 1 agcdn and the 
same proces repeats . Cyclic change in the modes of operation is 
given in Fig. 2. 2 . Ihe pattern of modes in steady state is 

1,2, 3, 2,1,.. . The typical wave-forms of Vci,Vcz,i8 in steady 
state are shown in Figs. 2.1(b). 


2.1.1 DIFFHiEMTIAL EQUATIONS OF TOE VOLTAGE DOUBLING KECTIFIHl 
FOR MODES 1,2,3 AND TOE METHOD OF SOLUTION 


In this section the differential equations describing the 



rectifier in different modes have been derived . The equations are 
solved for the expressions for voltage across capacitor 1, Vci(t), 
voltage across capacitor 2 , Vc 2 (t). The two variables of interest 
are Vci and Vc2 in this circuit . The reference directions for 
capacitor voltages and currents are shown in Fig. 2.1(a) . KVL and 
KCL have been applied to derive the differential equations . The 
diodes Di and D 2 are assumed to be ideal. 

NQRMRLISATIOH SCHEME .* The sinusoidal voltage source applied is 
Bnsinut. The arrplitude is taken as beise voltage . Base 

impedance is l/ojC . The base current thus becomes Bn.o.C. 
Normalised load ressistance as per this scheme becomes w.C.R. The 
equations derived for Vci,Vc 2 ,i 8 have been converted in p.u.form 
using above normalisation sdieme. This normalisation scheme helps 
to obtain design relations for the power doubling rectifier . All 
normalised variables are indicated by suffix n e.g. i«n is the 
normalised source current . 

Mode 1 (Fig. 2.3): In this mode Di is on and D 2 is off .Sotirce 
current is positive . The circuit enters Mode 1 at wt = a . The 
three differential equations are : 

source current, i# - C.EA.co.cos(ut) + ( Vci + Vc 2 )/R [2.1] 
capacitor 1 voltage ,Vci = E«.sin(ut) [2.2] 

(Vci + Vc2) / R = -C.d(Vc2)/dt [2.3] 

Equations 2.2 and 2.3 give 

( 1/R.C.+ D )Vc2 = -( En/RC ).sin(wt) [2.4] 

where D * d/dt 

Solution of above differential equation yields 

Vc 2 = -(Ei«)/( 1 + o^.R^.C^) .sin(wt) + Ki.exp( -t / RC ) 

+ (Bn.o).R.C)/(l + u^.R^.C^) .cos(ot) [2.5] 




Fig. 2. 2 Mode Transiticai diagram for the Voltage Doubling 

Rectifier fed by ideal sine wave voltage source 



Fig. 2. 3 Equivalent circuit in Mode 1 


Applying the initial conditions. 
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ut = a ,Vc2 = Vc 2 (a) in eq. 2.5, we have 

Ki = [ Vc2(a) + EVn.sin(a)/( 1 + w^.R^.C^)] .exp( a / w.R.C) 

- [ Bi,.Q.R.C.cos{a) / (1 4to^.R^.C^)].e3q)( o / w.R.C) 

[ 2 . 6 ] 

Normalised equations (as per normalisation scheme given above) for 
voltage across capacitor 1, voltage across capacitor 2 and source 
current ,i8 are given below : 

Vein = sin(wt) 

Vc2n = pl.sin(wt) + p2.cos(wt) + kln.eXp(-wt/Rn) 
isn - COS(wt) + (Vein + Vc2n)/Rn 

where Rn , pi , p2 and kin are given by 
Rn — W.R.C 
pi = -1/(R^ + 1) 

p2 = Rn/(I^ + 1) 

kin = kl/Bo 

Equations 2.4 , 2.5 and 2.6S 2.7 , 2.8 and 2.9 describe the 
operation of doubling rectifier in Mode 1. 

Mode 2 (Fig. 2.4) : This mode is characterized by the 

condition that both diodes are off and as a result source current 
is is zero . This mode starts at wt = p .The three equations are 
as follows : 


[2.7] 

[ 2 . 8 ] 
[2.9] 


source current , is 
capacitor 1 voltage , Vci 
capacitor 2 voltage , Vc 2 
where Vci(p) 


0 

Vc 1 0) . exp( -2/RC) ( t-p/w) 
Vc20)*.exp(-2/RC)(t-p/w) 
BB.sin(p) ( From Eq. 2.2) 


[ 2 . 10 ] 

[ 2 . 11 ] 

[ 2 . 12 ] 



Nonnalised equations are given below 
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Vein - sin(^).exp(-2/Rn)(ut-^) [2.13] 

Vc2n = Vc2(^)/B».exp(-2/Rn)(wt-^) [2.14] 

i.n = 0 [2.15] 

Equations 2.10, 2.11 and 2.12 & 2.13, 2.14 and 2.15 describe the 
doubling rectifier in Mode 2 . 

Mode 3 ( Fig. 2.5 } : iMs mode is diaracterized by the condition 
that diode Dz is on and diode Di is off .In other words source 
current ia is negative .Ibis mode starts at ut - s + a . The three 
equations for this mode are given below : 

capacitor 2 voltage ,Vc2 = -ElB.sin(wt) [2.16] 

source current , is » C.Eia.w.cosCwt) - ( Vd + Vc2 )/R [2.17] 

( Vci + Vc2 )/R = -C.d(Vci)/dt [2.18] 

Equations 2.16 and 2.18 give , 

Vci = EVi/(l + o^RV).[sin(«t)-«^RV.cos(ot)] + K 2 .exp(-t/HC) 

[2.19] 

* see i^pendix 1 for derivation of the expression for Vc 2 (a) 

Applying the initial conditions , wt = n + a , Vci - Vci(tc + a) 
in Eq. 2.19 , we have 
K2 = Vci(n + o)*+ 

B«/(l-Hi)^R^C^){sin(a)-ti)^R^C?.cos(a)}.exp( it + a )/(i)^R^C^ 

[ 2 . 20 ] 

Nonnalised equations for Voi , Vc2 and is are given below : 

Vein = l/(l+Rn^)[sin(ot)-Rn.COs((i)t)] + KZn .e^(-t»t/Rn) [2.21] 

Vc2n = -sin(«t) [2.22] 

isn = COS(Qt) - (Vein + Vc2n)/Rn [2.23] 

where lean = kz/Em 


Equations 2.16 , 2.17 and 2.18 & 2.21 , 2.22 and 2.23 describe 
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Method of soluticxi : The equations for different inodes derived 
above have been programned on HP 9000 using Fortran 77. The flow 
charts of the program are shown in Fig. 2. 6. The transient 
behaviour has been calculated from zero initial conditions and 
steady state haus been obtained after a no. of cycles. The initial 
conditions Vc2(a) ,Vc 2 (^) and Vci(n+«) have been calculated and 
given in Appendix 1. 


2.1.2.AVH»GE OOTPOT VOLmX 

In this section an expression for average output voltage 
has been derived .llie average output voltage is calculated over 
one half period i.e.n radians. Uie output voltage repeats every x 
radians i.e.the frequency of the ripple in the output voltage is 
twice the frequency of input supply . The expression for the 
output voltage has /thus ,been derived by integrating Voi 't Vcx 
from a to n t a in two intervals , oc to ^ and p to n * a . 


Average output voltage ,Vo = 1/n 


n*oc 

(Vc2 + Vci).d(wt) 


= 1/x 


a 

(Vci+Vc2)d(a)t)+ (Vci+Vc2)d(o>t) 

a 0 


[2.24] 

Substituting for Vci ,Vc2 from equations 2.2 , 2.5 , 2.11 and 2.12 
in 2.24 suitably , we get 


Vo = 1/jt { A + B + C + D) 


where A = B».[cos(a)-cos(p)] 

B = Rn/2.E4ii.8in(p).[l - exp{2p-2bt-2it)/w^R^C^ 


[2.25] 

[2.26] 
[2.27] 
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'Read Initial Conditions 
viz. Vci(0*), Vr.2(0*)and 
the value of the i 
Load Resistance / 


jt = 0 


Initialise time counter 


Is 

e ^ Vcl( 0 *) ? 


Ccnpute and store the values 
3 f Vci, Vc 2 , is and Vo(=Vct + Vc 2 ) 
fran eq. 2.9, 2.12, 2.13 
with current initial conditions 


(i)t - ojt + tstep 


/ Is ' 

ut ^ jr 


s® ^ Vcl ? 


Is \ NO Is N 

e ^ Vc 2 = ? >>1 V^t = 2 ji 


npute and store the values 
of Vci, Vc 2 , is and Vo 
rom eq. 2.6, 2.7 and 2.8 
ith current initial conditions 


Vo(n+l)rp - Vo(„)Tp = £. 

\ £ = 0.01 


(i)t - cat + tstep 


Corpute and store the values 
of Vci, Vc2, is and Vo 
from eq. 2.20, 2.21 and 2.22 
with current initial condition; 


A COKPUie AM£> 

/ PRINT AvtRAGiH 

'and f-To-?K\ftl£ 


ut - (jt + tstep 


tl. 






C = pl.Ei».[cos(a)-cosO)]+p2.Eiii. [sin(p)] 

- p2.Ei«.[sin(a)]-Ki .w^RV.[exp(-^/uVc^)-exp(-a/wVc^)] 

[2.28] 

D = uVc^/2.Vc20).[l-exp(2{a+7r}-2^)/wVc^] [2.29] 

Von - Vo / B* * (1/«).B,.( A + B + C + D ) 

Ion - Von/Rn 

Derivations of expressions for A,B,C and D are given in A^^^endix 1 

2.1.3. PEAK TO PEAK ODTTOT VOLTAGE RIPPLE 

In this section an expression for peak to peak output 
voltage ripple has been derived. The output voltage follows a 
certain pattern .It reaches maxiirun at ut s p and touches the 
lowest point at ut a: a. This is the period when diode Di is on and 
D 2 is off . The same voltage rii^le is observed during the period 
from ot = jt + a to ot = K + p, the period during vMch diode Dz is 
on Di is off .Hence peak to peak voltage ripple is calculated 
during this period only. 

Vrippie(p-p) = - (Vo)^^^ 

Vripple(p-p) = ~ 

Equations 2.2 and 2.5 in 2.31 give 

Vripple(p-p) = EW.sinO) + [Ei,.W.R.C./(l4toVc^).COsO)] + 
[Kl.exp(-^/w.R.C)]- E«/(Hw^R^cf).sinO) - B«.sin(a) - Vc 2 (a) 

[2. 31a] 

Normalisation gives 

Vripple(p-p)n — sin{/5)+ [Rn/(1+Rn^) .COS{^)] + [Kin .e35*(~P/R*»)] 

- l/(l+Rn^).sinO) - sin(o) - Vc2n(o) [2.3Ib] 



The expressicxi for peak to peak voltage ripple in original and 
normalised form could be used for design purposes. 

2. 1.4. DISCUSSION OF SIMULATION RESULTS 

In this section the general behaviour of the doubling 
rectifier has been studied . The transient behaviour, steady state 
bdiavi our, average output voltage and peak to peak voltage ripple 
for different values of load resistances have been calculated and 
shown in Figs. 2. 7, 2.8 and 2. 9. These results are discussed below : 


2. 1.4.1 Loaded Trai>sient and Steady State B^iaviour 

Pig. 2.7 shows the sinulated output voltage, volt age across 
capacitor 1, voltage across capacitor 2, source current is and 
input voltage for the doubling rectifier with load resistance of 
10.0 p.u. 

The circuit of figure 2.1 is excited fron cold and its 
transient behaviour is noticed . At «t=0^ , diode Di turns on as 
the supply voltage is larger than voltage across capacitor 1 
(being Ol.This allows the capacitor to follow the supply voltage 
(the circuit is in Mode 1 fdiere Vein = sin(ot)) .The voltage 
across capacitor 2 increases in the negative direction(see Fig 
2.3) .The circuit renains in Mode 1 till source current is greater 
than zero .The monnent it crosses zero and tries to increase in the 
negrative direction ,the diode Di ceases conduction and the circuit 
transits to Mode 2 idiere in both the diodes are off (Fig. 2.4). 

Now the capacitor 1 discharges and charges capacitor 2 
slowly tdiich was hitherto holding negative voltage across its 
terminals . The condition of both diodes off remain till a point 
idien the voltage across capacitor 2 becomes equal to negative of 
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voltage across OUTP.UT,CAPl.CAPZ (p.u.) Rrjo o p.v^. 



OM <ao om OM wt QM lua oia QOS OM OM 


TlMlE(sec) 



TIHE(sec) 


Fig. 2.7 Loaded Transient Behaviour of the voltage doubling 
rectifier fed by ideal sine wave voltage source 



the siiM>ly voltage . At this point Dj_ turns on renains off 
(since voltage across it is greater than supply voltage ) and the 
circuit transits to Mode 3 .The output voltage , capacitor 2 
voltage increase and capacitor 1 voltage decreases. 

ligain the circuit transits to Mode 2 vrtien the source 
current crosses zero . Prom here on the diodes Di and Dz remain 
off till a point tdten voltage across capacitor 1 becomes less than 
st^ly voltage . 

The circuit goes through above steps in transient and 
afterviards the voltage across capacitor 1 , capacitor 2 aiid output 
follow a periodic pattern and the circuit moves in a cyclic order 
from Mode 1 to Mode 2 to Mode 3 to Mode 2 to Mode 1 and so on . 

2. 1.4. 2 Average Output Voltage Vs. Load Resistance (p.u. ) : 

The characteristics is shown in Pig. 2.8 .Averaige output voltage 
rises as the load resistance increases . The output voltage 
reaches a value of approximately 2.0 p.u. at substantially high 
values of resistances. With the increased load resistances , the 
capacitors once charged to 1.0 p.u. do not discharge much in Di 
and D 2 off mode and as a result output average voltage goes vip . 
It is seen from Fig. 2.8 that an acc^table voltage doubler with 
average output voltage greater than 1.9 p.u. must have a load 
resistance greater than approximately 100.0 p.u. Thus value of 
capacitance can be calculated for a given load and frequency. 

2. 1.4. 3 Peak to Peak Output Voltage Ripple Vs. Load Resistance ; 

The characteristics is shown in Pig. 2.8 . The peak to peak 

voltage ripple reduces as we go on increasing the load resistance. 
The ripple reducss owing to the fact that capacitors do not Icse 
their charge heavily as conpared to the case when load resistance 
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aVERRGE OUTPUT VOLTAGE vs LORD RESISTRNCE(p^ 
FOR SQUARE WAVE AND SINE WAVE CASE 



Average CXitput Voltage Vs Load Resistance 
characterstics for the voltage doubling rectifier fed 
by ideal sine wave voltage source 


Fig.2.8 



PERK TO PERK Va.TRGE RIPPLE FOR 
SQURRE mVE RND SI^C mVE VOLTRGE SOURCE 
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Fig. 2. 9 


Peak to Peak output voltage Ripple Vs Load Resistance 
characters tics of the voltage doijbling rectifier fed 
by ideal sine wave voltage source 



15.8 p.u. 
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2. 1.4. 4 Average output current follows the same pattern as that 
of average output voltage for average output voltage is the 
product of average output current and load resistance . 

2. 1.4. 5 When the load resistance is theoretically increased to 
infinity i.e. the output terminals of the power doubling rectifier 
are open circuited ,the following happens . Assianing we started 
with zero initial conditions ,in the first half cycle till ot = 
n/2 ,the capacitor 1 voltage will follow the ii^t voltage and 
capacitor 2 voltage remains at zero voltage . At ot = Tt/2 , the 
diode Di goes off and capacitor 1 remains at 1.0 p.u. for all 
values of time after it as there is no path to discharge . At wt = 
n , Dz turns on voltage acaross cjapacitor 2 will be 180** out of 
phase with input voltage . At ot = 3it/2, Dz goes off and capacitor 
2 stays at 1.0 p.u. .The output voltage becomes 2.0 p.u. and 
remains at this valtie for any further values of time . 

2.2 VOLTAGE D0DBLIN8 RBOTIFIER FED BY HON Ite2U:i SINE NAVE VOLTAGE 

SOURCE 

The doubling rechifier ( Fig. 2.10 ) fed by nonideal sine wave 
voltaige source has three inodes of operation . The non ideal sine 
wave voltage source has resistance , R« and reactance, uL« . The 
inodes are characterized by patterns of diode cxnduchion . The 
modes of cjperation are 
. . Mode 1 :Di ON , Dz OFF, a< ut i 

Mode 2 :Di OFF, Dz OFF, ^<<^tSJt + aandIt + ^<uts2rc+a 
Mode 3 :Di OFF, DzON,n+a<ot:Sir + p 




except for the fact that the voltage across capacitor 1 is not 25 

same as the svqply voltage in mode 1 and voltage across capacitor 2 
is not same as the negative of the supply voltage in mode 3. This 
is so because there is some voltage drop across source irtpendence. 

2.2.1 DIFFEROrriAL BQCnTIONS FOR MODES 1, 2 AND 3 AND HIE METHOD 
OF SOUJTICN 

In this section the differential equations describing the 
rectifier in different inodes have been written . The derivatives 
of voltage across capacitors 1 S 2 and source current have been 
used in simulating the b^iaviour of the doubling rectifier . Hie 
three variables of interest are Vci , Vc 2 and i« . The reference 
direction of capacitor voltage and current and source current are 
shown in the Fig. 2.10 . The equations have been normalised as per 
the normalisation scheme given in the beginning of the chapter . 

The diodes have been assunned to be ideal. The pattern of modes in 
steady state is 1,2, 3, 2, 1,2... . The mode transition 

diagram is given in Fig. 2.11 . 

Mode 1 (Fig. 2.12): In this mode Di is on and D 2 is off 
current is positive . The circuit enters Mode 1 at ut = 
three differential ec[uations are 

E«.sin(<i>t) = i«.R8 + L«.d(i«)/dt + Vci 
(Vci + Vcz) / R = “C.d(Vc2)/dt 
source current, i* * C.d(Vci)/dt + ( Vci + Vc2 )/R 
Normalising the above equations , we have 
(from equation 2,32) 

sin((ot) = U.IU/Em + (L8/B«).[d(is)/dt] + Vci/Ei» 
sin(«t) = (is.Rs .(■)C)/(E«.(i)C) + (L«.<i)^C)/(Ei«.{i)C).[d(i8)/d(«t)] 


. Source 
a . The 

[2.32] 

[2.33] 

[2.34] 








sin(<i)t) = isn.Rsn + ((i)L« )n . 


+ Vein 


[2.35] 


(frem equaticxi 2.33) 

( Vcl + Vc2 )/Eb = - (RC/En).d(Vc2)/dt 

Vein + Ve2n = - (RC.0))/(&.).d(Ve2)/d(wt) 
Vo,n+Vc2„ = -R„.[-|^] 

{from equation 2.34) 


[2.36] 


i«/(B«.oC) = C/(&i.a)C).d(Vci)/dt + (Vci + Vc2)/(R.E,.oC) 


isn m 


d(Vcl )nl 

■3(0)1) J 


(Vein + Vc2n) 


[2.37] 


In all the eqruations derived above suffix n r^resents normalised 
variables e.g. (o>Ls)n = o) .L«.C represents normalised source 
reactance . 

Equations 2.32 , 2.33 , 2.34 , 2.35 , 2.36 and 2.37 describe the 
operation of voltage doubling rectifier in Mode 1 . 


MCX£ 2 (Fig. 2.13) : In this mode both the diodes Di and Dz are 
off . Hie circuit is characterised by source current , is = 0 .Hie 
circuit enters mode 2 at oit = p tdien source current becomes zero . 


Hie three differential equations are 

source current /is = 0 [2.38] 

(Vcl + Vc2) / R = -C.d(Vc2)/dt [2.39] 

(Vcl + Vc2) / R = -C.d(Vci)/dt [2.40] 

Normalised equations can be written as follows 

isn “ 0 [2.41] 

Equations 2.39 and 2.40 are of the kind of equation 2.33 , hence 
the normalised equations directly become 

Vein + Vc2n - “ Rn . ^ d(ot ) j [2.42] 
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Egxiations 2.38 , 2.39 , 2.40 , 2.41 , 2.42 and 2.43 describe the 
operaticm of the dcfubling rectifier in mode 2 . 

MCIS 3 (Fig. 2.14) : In this mode diode Di is off and diode D 2 is 
off . The circuit alters mode 3 at wt = n + a »dien negative of the 
supply voltage beccjmes greater than voltage across capacitor 2 . 

The source current is negative . The three differential equations 
of the circuit are derived below 

Ek.sin(wt) = i..R* + L«.d(i8)/dt - Vc2 [2.44] 

(Vci + Vc2) / R = - C.d(Vci)/dt [2.45] 

source current, i« = - C.d(Vc2)/dt - ( Vci + Vc2 )/R [2.46] 



Equaticxis 2.44 , 2.45 , 2.46 , 2.47 , 2.48 and 2.49 describe the 
operation of the doubling rectifier in mode 3 . 

MEmOD OF SCXjUTION : The equations derived above for the 
operation of the doubling rectifier in different modes have been 
programmed on HP-OX 9000 in Fortran 77 . Hie expressions for the 
first derivatives of the variables Vci , Vc2 and i« are required . 
Hie three simultaneous equations in variables Vci , Vc2 and is are 



Appendix 3 gives the details of this subroutine. 

In mode 1 the corrputed values of the source current are 
continuously mcxiitored aixi as soon as the zero crossing of source 
current is detected , the program flow is either transferred to 
corputations in Mode 2 or Mode 3 . Normally the program flow is 
transferred to corputations in Mode 2 . 

In mode 2 the corputed values of the variables Vci and 
Vc 2 are continuosly monitored .Ihe program flow is either 
transferred to Mode 1 or Mode 3 depending on whether the supply 
voltage is greater than voltage across capacitor 1 or the negative 
of the sijpply voltage is greater than voltage across capacitor 2 . 

In mode 3 the conputed values of the source current are 
continuously monitored and as soon as the zero crossing of source 
current is detected , the program flow is either transferred to 
corputations in Mode 2 or Mode 3 . Normally the program flow is 
transferred to conputations in Mode 2 . 

fihe conditions which decide these transfers are already 
described in the beginning of this chapter .Flow chart of Fig. 2. 15 
further describes this scheme .'Hie conputed values of the 
variables Vci , Vc2 and is are stored appropriately to be used 
later for conputation of average output voltage and peak to peak 
output voltage ripple . 

2.2.2 MEIHOO OF OGMHJTATION OF AVERAGE OUTPUT VOLTAGE AND PEAK TO 

PEAK OUTPUT VOLTAGE RIPPLE 

This is conputed in the manner as follov^ : 

The data points are available every 'tstep' radians . 'tstep' is 
the step size in conputation . We assume the output voltage to 
remain constant over 'tstep* radians .The smaller the value of 
'tstep* , higher will be the computational accuracy . Hence tte 
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t = 0 


Vcl(0*) 


fis. 3.i5<<») 

confutations according to 
the Flow Chart for Mode 1 


confutations according to 
the Flow Chart for Mode 3 

^l<l 2-15tC> 


fl&.ZtSCb) 

confutations according to 
the Flow Chart for Mode 'i 


go to mode 1, mode 2 or mode 3 
depending upon the value of 
source current, is 


Using stored values of 

confute and Print 

numerical Average output voltage and 
Peak to Peak voltage Ripple 


ia and Vo 


Flow chart for the Similation of the Voltage 
doubling rectifier fed by nm ideal aine wave source 













store the values 

Cornpute ^ Vo(=Vo. + Vc2) 

of V.,, V.2, ^ ,.3,C3.^7j 

from eq. 2.35^*^3re wave} voltage sor 
for sine wave(sq _. conditions 

with current i 


,orce 






t- LOW Ui fl^K I KUK r^Oi^v^ z 


Compute and store the values 
of Vc.i, Vc/., is and Vo(=Vci + Veil) 

frail eq. 2. m(3-3j),2.l42.(3.3^ , 2-45 (3.3 3j 

for sine wave(square wavej voltage sorce 
with current initial conditions 



ut = (i)t + tstep 


Is 

ut ^ 2 jc ?, 


end of cycle 


Is 

e Vc2 = 


Is \ 
e ^ Val = ? 



Modes 



1 change initial conditions 
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V / 


do TO 
MODE / 1 

<so To 
MOPE 3 
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PLOU) CH/V^T Po^ KODB" 5 



Go TO 

HODEa 


e ^ Vci 


<5o TO 
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average output voltage in per unit is 


\T _ V* (Vclnl + Vc2nl ) 

VoBvon - 2^ i— 

1«1 “ 

where k = no. of data points over a period of 2n radians. 
In the coniplete process of transient and steady state there can be 
p.k data points vfliere p is the no. of cycles taken by the circuit 
in reaching the steady state. Only last k datapoints, however, 
are retained for these are the data points corresponding to the 
steady state of the circuit. The data points during the transient 
state are only needed for plotting . 

The naxinum and imniirun values of the output voltage are picked 
and the difference is outputed as peak to peak output voltage 
ripple . 

2.2.3 DISCDSSICM OF THE SDfJIATIOK RESOLTS 

In this section the general bdiaviour of the doubling 
rectifier las been studied . Hie transient brfiavi our, steady state 
behaviour, average output voltage and peak to peak voltage ripple 
for different values of load resistances and source reactances 
have been calculated and shown in Figs. 2. 16, 2.17 aixl 2.18. 

These results are discussed below : 

2. 2. 3.1 Loaded transient ax>d stead^T state b^viour of the 
doubling rectifier 

•Hie brfjaviour is shown in Fig. 2.16 . The doubling rectifier 
with load resistance of 10.0 p.u. , source resistance of 0.01 p.u. 
and source inductive reactance of 0.05 p.u. was simulated for 
studying its transient and steady state behaviour . 

Fig. 2.16 shows the simulated output voltage , voltage 
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Fig. 2. 16 


Lcsaded Transient Bebaviour for the voltage doubling 
rectifier fed by acn ideal sine vrave voltage source 







and input voltage .The doubling rectifier is excited from cold and 
its transient behaviour is recorded .At ut = 0^ , diode Di turns on 
as the supply voltage is larger than the voltage across capacitor 
1 (being zero). Source current, i* jumps to a value decided by the 
series combination of R-L-C. Source current is a sinusoidal pulse. 
The diode Da is off as it is reverse bicsed. The voltage across 
capacitor 2 increases in the negative direction. The circuit 
remains in Mode 1 till source current is greater than zero .Hie 
moment it crosses zero and tries to increase in the negative 
direction, the diode Di ceauses conduction and the circuit transits 
to Mode 2 vdiere in both the diodes are off. 

Now the capacitor 1 discharges and charges capacitor 2 
slowly vhich was hitherto holding negative voltage across its 
terminals. The ccxidition of both diodes off renain till a point 
when the voltage across capacitor 2 becomes equal to negative of 
the supply voltage. At this point Dz turns on Di remains off 
(since voltage across it is greater thaoi si^ly voltage ) and the 
circuit transits to Mode 3. Hie output voltage, capacitor 2 
voltage increase and capacitor 1 voltage decreases. 

Again the circuit transits to Mode 2 idien the source 
current becomes zero. Prom here on the diodes Di and Dz remain off 
till a point when voltage across capacitor 1 becomes less than 
stpply voltage . 

Hie circuit goes through above steps in transient and 
afterwards the voltage across capacitor 1 , capacitor 2 and output 
follow a periodic pattern and the circuit moves in a cyclic order 
from Mode 1 to Mode 2 to Mode 3 to Mode 2 to Mode 1 and so on. The 
mode transition diagram is shown in Fig. 2. 11. 


2. 2. 3. 4 PERIODICITY OF THE COTPUT Va/HYIE HAVEFOSM 



The siirulation results were recorded for various of load 
resistances and a very peculiar behaviour of the doubling 
rectifier was noticed. 

Normally the output voltage has a ripple frequency twice ais 
that of the input supply frequency. But, when the doubling 
rectifier is fed by a nonideal sinusoidal voltage source, the 
output voltage has a ripple frequency twice as that of the input 
supply voltage frequency till a critical value of the load 
resistance. Hie critical value is 30.61 p.u. when the source 
reactance is 0.05 p.u., source resistance is 0.01 p.u. and load 
being supplied is 10.0 p.u. Hie frequency of the output voltage 
keeps on reducing as the value of load resistance is increased 
beyond 30.61 p.u. At very large values of load resistances the 
ripple frequency is so low that for all practical purposes the 
output voltage can be assimed to be flat. Hie table 2.1 shows the 
growth in the time period of the ripple in the output voltage. 

2. 2. 3. 2 VARIATION OP TBE AVERAGE OOTFOT VOLTAGE WIHI LORD 

RESISTANCE AND SOORCE REACTANCE 

Average output voltage vs. source reactance : 

Hie characteristics is shown in Pig. 2.17. From Fig. it is clear 
that the average output voltage falls both at very low values of 
uLa and at at relatively large values of (oLs . Hie average output 
voltage is close to the ideal value of 2.0 p.u. for the range of 
oiLa ~ 0.015 to 0.07. 

From the wave forms of the Fig. 2. 16, the source current is 
seen to have sharp sinusoidal pulses at uLs = 0.05 p.u. As wLa 
reduces towards ideal value of zero, the araplitixle of the pulses 
become very large as analysed in section in 2.1. Due to th^e 
sharp pulses of source current, the diodes Di and Dz turn off 
affoi- rvunAintino for very short intervals of time in positive and 



TABLE 2.1 

VARIATICW OF RIPPLE PREQqENCy OF THE OJTPUT VOLTACTl WITH 
VARIATICW IN LORD RESISTANCE 


Load Resistance 

(p.u.) 

Time period as a multiple of 

of input supply time period 

i. 

10.4) 

0.5 

2. 

25.0 

0.5 

3. 

30.61 

0.5 

4. 

30.62 

a.o 

5. 

39.32 

2.0 

6. 

39.33 

4.0 

7. 

100.51 

11.0 

8. 

100.52 

26.0 

9. 

1000.0 

very large(see fig. 2. 20) 





Pig. 2. 17 ’’Average Output Voltage Vs Source Reactance for 

varying Load Resistances characterstics" for the 
voltage do\jbling rectifier fed by non ideal sine 
vrave voltage source 



(*n* d) 36 b!H 0A :jnd:5.no 
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Fig. 2. 20 Loaded Transient behaviour for the voltage doubling 

rectifier fed by ncai ideal sine wave voltage source 
at a load resistance of 1000.0 p.u. 
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negative half cycles respectively. Therefore the capacitors sv^jply 
energy to the load for most part of the cycle and thus output 
voltage falls for most part of the cycle. The average output 
voltage, hence, is low at low values of wLs. 

At high values of wLs ( > 0.07 p.u.) the peak sotirce current 
reduces and due to large values of the source reactance the 
conduction period of the diodes in the positive and negative half 
cycles Increases. However, the average output voltage falls as the 
inductive circuit is not able to transfer charge quickly. 

Average output voltage vs. load resistance : 

The b^iaviour of the output average voltage vs. load resistance is 
quite complex in nature. The average output voltage rises as the 
load resistance is increased for the values of o)L« in the range 

from 0.015 to 0.07. At higher valvies of the source reactances 
(>0.07 p.u.), the average output voltage rises till a critical 
value of the load resistance. Increase in the load resistances 
beyond this value causes a decrease in the average output voltage. 
For small values of the source reactance, as the load resistance 
is increased, the capacitors once charged do not lose much of 
their stored energy and as a result the output voltage rises. 
However, at large values of the source reactances ( >0.07 p.u.), 

the increase in the load resistance causes a decrease in the 
average output voltage becaijse the energy is trapped in the source 
inductance. 

2. 2. 3. 3 VARIATION OF PSUC TO PEMC OUTEOT VOLTAGE RIPPLE WITH MfiD 


RESISTANCE AND SOORCE REACTANCE 






Peak to Peak Output Voltage ripple vs. Source reactance: 
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The cliaracteristics is shown in Fig. 2.18. Jls explained in 
sec. 2. 2. 3. 2 above, the ripple voltage increaises as the source 
reactance is increaused. At low values of the uL*, the diodes 
conduct large current pulses for shorter durations. These large 
current pulses cause sharp voltage rises across capacitors IS 2. 
For most part of the cycle diodes are off and capacitors discharge 
from their peak levels causing a large peak to peak output voltage 
ripple. 

At large values of source reactances, the conduction period 
of the diodes increase and due to the continuous change in the 
supply voltage ( the input sinusoidal voltage magnitvide is 
continuously changing ) the ripple in output voltage increases. 
The bdiaviour of the percentage ripple is conplex due to changes 
in the average output voltage with source reactance(fig. 2.17). 

Peak to Peak output voltage ripple vs. load resistance: 

The behaviour of the peak to peak output voltage ripple is inverse 
of the b^iaviour of the average output voltage, as seen in fig. 


2.18. 



CHAPTER 3 


VOLTAtffi DOOBLING IfflCTIFIER FBD BY SQUARE WAVE VM.TAGE SOURCE 

In this chapter analysis and siniulatioti of the voltage 
doubling rectifier fed by ideal and non ideal square wave voltage 
source has been done. Simulation results have been discussed. Some 
iinportant conclusions have been made. 

3.1 VC»:.TA(3E DOUBLING RBCTIFIHl FH> BY 

AN IDEAL SQUARE WAVE VOLI^GB SOURCE 
Ihe doubling rectifier of Fig. 3.1 with ideal square wave 
voltage source has two inodes of operation . These inodes 2 ure 
characterised by the patterns of diode conduction . The mode 
transitions are shown in Fig. 3.2 .The inodes are : 

Mode 1 : Di on , D 2 off, 0 < wt < n. 

Mode 2 : Di off , Da on , ji < ut < 2ic. 

The description of each mode is as follows : 

1 : The equivalent circuit for this mode is shown in Fig. 3. 3. 
In this mode the diode Di starts conduction at ot = 0 as the 
svqpply voltage is larger than capacitor 1 voltage . At ut = 0 ,an 
inpulse of current flows through diode Di as the voltage across 
capacitor 1 is less than E« .Voltage across capacitor now juips to 







EVk from a value less than Bn as a result of this infnilse of 
current .Had there not been any parallel path across capacitor 1, 
diode Di would have stopped conduction for the source current 
would then have become zero .Since we have a series R-C in 
parallel with capacitor 1 , the diode Di keeps on conducting to 
supply charge to capacitor 2 and hence the source current ( same 
as diode Dt current) is not zero .Voltage across capacitor Ca 
decreauses in Mode 1 . At wt = jt , the sipply voltage jutps to -Bn 
and the circuit transits to Mode 2 due to the following : 

At o)t = It" ,the diode Di is on and diode Dz is off . At wt = it, 
the source voltage jutps to -Bn .'Hie voltage across diode Di 
reverses and becomes -2Ein .This causes the source current to 
reverse aiKl as a result ttim diode Di off . Diode Da is forward 
biaused as the negative of the suqpply voltage is laurger than the 
voltage across capacitor 2 . This turns diode Da on . 

MM® 2 ; The equivalent circuit for this mode is shown in Pig. 3. 4. 
At at = It* , the diode Da ttims on as the negative of the supply 
voltage is larger than the voltage across capacitor 2 . Source 
experiences a negative inpulse of current at wt = n* . This is due 
the difference(AV) in voltage across capacitor 2 and negative of 
supply voltage at wt = it*. The strength of the inpulse is 

(C.AV). Larger the difference larger will be the strength of the 
inpulse. The diode Da remains on till «t = 2it tdien Di turns on 
oning to the reasons ejplained above in Mode 1 .The circuit 
transits to Mode 1 at wt = 2it .Cyclic change in modes of operation 
is shown in Fig. 3.2 . The patterns of modes in steady state is 



1, 2, 1, 2.. . The typical waveforms of Vci, Vc 2 and i* are given 
in Fig. 3.1(b) . 

3.1.1 DIPFTRENTIAL BQOATIONS OF TEffi DOUBLIMO RECTIPIBR FOR MODES 
1 m> 2 m> TBE MSmOD of SOLDTIOM 

In this section the differential equations describing the 
rectifier in different modes have been derived . The equations 
are solved for the expressions of Vci and Vc 2 . The two variables 
of interest are Vci and Vc2 in this circuit . The reference 
directions for capacitor voltages and currents are shown in 
Fig. 3.1 . KVL and KCL have been applied to derive differential 
equatioaos . The diodes Di and D 2 have ben assumed to be ideal . 

Normal isaticxi Scheme : Normalisation scheme is same as that 
given in chapter 2 (Sec 2.1.1) .The equations have been converted 
in p.u. form using this normalisation scheme .All normalised 
variables have been shown with subscripts n . 

MODE 1 (Fig. 3.3) : In this mode Di is on and D 2 is off . Source 
current is positive . The three differential equations in Mode 1 
are given below. The circuit enters Mode 1 at »t = 0*. 
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Fig. 3. 2 Mode Transition Diagram for the voltage doubling 

rectifier fed by an ideal square wave voltage source 




Fig. 3. 4 


Equivalent Circuit in Mode 2. 
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Fig. 3. 2 Mode Transition Diagram for the voltage doubling 

rectifier fed by an ideal square wave voltage sovirce 



Fig. 3. 4 Eqpii valent CLrcmt in Mode 2. 





Vcl = Em 
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( Vcl + Vc2)/R = -C.d(Vc2)/dt 
Sotirce current 4* = (Vci + Vc2)/R 
[3.1] and [3.2] give 

( 1/RC + D )Vc2 = -Em/EC ; where D ■ d/dt 

Solution of above differential equation gives 

Vc2 = -& + [(Vc2(0*) + Em ].exp(-o)t/o.RC) [3.4] 

vdiere Vc 2 { 0 *) is the voltage across ceq^citor 2 at wt = 0* . 
Normalised equations as per normlisation scheme discussed above 
for Vcl ,Vc2 and i» zure given below : 


[3.1] 

[3.2] 

[3.3] 


Vein = 1.0 [3.5] 

Vc2n = -1 + [Vc2n(0'*') + l].e^(-ti)t/Rn) [3.6] 

where Rn = <aRC (normalised load resistance) 

isn — (Vein + Vc2n)/Rn [3.7] 

In steady state the value of Vc2(0*) is . Hence equation 
3.6 changes to the following 

Vc2n = -1 + 2.e3q?(-<i)t/Rn) [3.8] 


Equations 3.1 , 3.4 , 3.5 and 3.6 describe the operation of 

voltage doubling rectifier in Mode 1 . 

MXS 2 (Fig. 3.4) : This mode is characterised by the condition 
that diodes Di is off and D 2 is on . In other words source current 
,i» is negative .This mode starts at ot = n* .The three 
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differential equations under this condition are given below : 

Vc2 - Em [3.9] 

( Vci + Vc2 )/R = -C.d(Vci)/dt [3.10] 

i* = -(Vci + Vc2)/R [3.11] 

Equations 3.9 aiKl 3.10 give 

( D + 1/RC )Vci = - B./RC ; where D ■ d/dt 

Inlying initial condititons that Vci = Vci(n) at «t = K*and 
solving ,we have 


Vci =-£. + [&+ Vci(jt*)].e3ip {-(0t-it)/<dRC} [3.12] 

In steauiy state Vci(jt*) = Bin .Hence equation. 3.12 modifies to 

Vci = -Em + 2.EV..exp {-(«t-ji)/uRC} [3.13] 

Normalised equations for Vci and Vc2 are griven below 

Vein = -1 + [ 1 + Vcln(jl*) ].exp {-(ut-n)/Rn} [3.14] 

Vc2n = 1.0 [3.15] 

i»n = -(Vein + Vc2n)/Rn [3.16] 

In steady state the eq.3.14 reduces the following 
Vein = -1 + 2. exp {-((i)t-Jt)/Rn} 


Equations 3.9 , 3.11 , 3.13 ,3.14 , 3.15 aiid 3.16 describe the 
operation doubling rectifier in Mode 2 . 

Method of Solution : The eqiations for different modes derived 
above have been programmed on HP-OX 9000 using Fortran 77 

.:eimtral library 

i I. T., KANPUR 



language. Hie flow chart of this program is shown in Fig. 3.5. 
The transients Izist for one half cycle only and steady state is 
reached iimediately after that due to the ideal natiure of the 
source . 

3.1.2 AVERAGE OOTFOT VCA.T7V3E In this section ,an expression for 
average output voltage has been derived . The eiqpression for 
average output voltage has been derived by integrating (Vci + 
Vc 2 ) over a period from 0 to 2ir . 

Average output voltage , 

zn 

Vo = l/2ft.J(Vci + Vc2).d(wt) 

0 

It 2lt 

=l/2rc.[ J(Vci + Vc2).d(ut)+ J(Vci +Vc2).d(ot)] 

0 K 

Substitution from equations 3.1 ,3.4 ,3.9 ,3.13 gives in above 
equation 

Vo = (l/n).2B«.0).RC [ 1 -e3qp(-Jt/(oRC)] [3.17] 

Normalisation gives 

Von = 2Rn/rc.[ 1 - exp(-Jt/Rn)] [3.18] 

Average output current , 

Ion = 2/ir.[ 1 - e3qp(-Ji/Rn)][3.19] 

Equations 3.18 and 3.19 give average output voltage and current . 

3.1.3 PEAK TO PEAK OOTTOT VOLTNX RIPPLE In this section an 
expression has been derived for peak to peak output voltage 
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ripple. Hie output voltage follows a certainpattem .It reaches 
inaximjfn at ut = Ot n* 2n* . . and touches the lowest point at wt = 
0 , n , 2n. . . Hiat is why output voltage ripple can be calculated 
in the period from 0 to n or ji to 2rc . Typical output voltage 
waveform has been shown in Pig. 3. 1(b). 

Hence Peak to Peak output voltage ripple, 

Vo-ripple pp = (Vo) ♦ - (Vo) - 


= (Vc2) + - (Vc2) 

‘^ut -0 «t«n 

From equation 3.4 ,we have 

= 2E1»[ 1 - e3^(-jt/(i)RC) ] 

Normalised peak to peak voltage ripple, 

Vo-ripple ppn = 2.[ 1 ” e3!p(”1t/Rn) ] 


[3.20] 


[3.21] 


With all the expressions found above it is clear that variables of 
interest in the power doubling rectifier , average output voltage 
and peak to peak output voltage ripple depend t{>on the value of 
normalised load resistance . 


3.1.4 SCONCE CURRENT IMFOLSE In this section , an expression for 
source current impulse has been found .At ot = 0,n,2x .. there are 
positive and negative source current impulses in order to svq^ly 
charge to capacitor 1 S 2 respectively . These impulses raise the 
voltage across capacitor 1 & 2 to Ein .The strength of impulse 
could be found as follows : 

Impulse strength = C.AVcin 


= C.AVc2n 



From the expession for peak to peak output voltage ripple found 
above, Iinpulse strength = 2.[ 1 - exp(-n/Rn) ] [3.22] 


3.1.5 DISCXJSSION Of SmjLATICXI RESULTS In this section the general 
behaviour of the rectifier has been studied .Ihe transient 
behaviour , steady state brfiaviour , average output voltage and 
peak to peak output voltage ripple has been recorded for different 
values of load resistances .Fig. 3. 6, 2.8 and 2.9 show these 
respectively. The results are discussed below : 

3. 1.5.1 Loaded Transient and Steady State B^iaviour 

Fig. 3. 6 shows the simulated output voltage , voltage 
across capacitor 1 and 2 and source current ,i« for the doubling 
rectifier with load resistance of 10.0 p.u. 

The circuit of Fig. 3.1 is excited from cold its b^iaviour 
recorded . At wt = 0*,the diode Di turns on as the sxipply voltage 
(+E1») is larger than capacitor 1 voltage (being zero). A source 
current impulse results because of this voltage difference .This 
voltage difference is [ - Vci«i)t)«o” ].The diode D 2 is off as 
the negative of supply voltage is less than voltage across 
capacitor 2 . With Di on and D 2 off the following happens : 

The voltage across capacitor 1 stays at and voltage across 
capacitor 2 starts falling from 0 to a value decided by RC .At wt 
= jt ,the source voltage jtiips to -Eii and diode Di tiims off . The 
diode D 2 turns on as the negative of the svpply voltage is greater 
than voltage across capacitor 2 .A negative source current inpulse 



sou'ce cun*ont{pjj.) ^ capi.cap?,cx.'t.pi-*w(p.uWaili'a;ge'»- 


R(p.a)=10.0 



Fig. 3. 6 


Loaded Transient Beh4\^our of 
rectifier fed by an ide?rl square 


•rxa 

► 



the voltage doubling 
rave voltage source 


flows throxigh D 2 to sv 5 >ply the charge required by the capacitor 2 
to acquire the voltage of En .A sindlar positive source current 
inpulse resulted at ut = 0* to supply the charge to capacitor 1 to 
enable it acquire the voltage of +E« .Hie diode Di remains off 
during the period from n to 2n .Hie diode Dz remains off during 
the period f ran 0 to n . 

At ut = 2ji , the diode Dz turns off as the supply voltage 
junps to -i-Em from -Bn and from here onwards the same process 
repeats . So the circuit follows a sequence of Di on ,Dz off and 
Dz on ,Di off and so on .Mode transition diagram is shown in 
Fig. 3. 2 . 

3. 1.5. 2 Average Output Voltage vs. Load Resistance (p.u.) 

Average output voltage rises as the load resistance 

increases .Hie output willreadi a itaxiirum of 2.0 p.u. with load 
resistance infinity , theoretically .With the increased load 
resistance , the capacitor 1 S 2 once diarged to 1.0 p.u. at ut = 
0*& n* respectively do not discharge irudi and as a result average 
output voltage goes rp .Fig. 2. 8 shows this bdiaviour . 

3. 1.5. 3 Peak to Peak Output Voltage Ripple vs. Load Resistance 
Hie peak to peak output voltage ripple decreases as the 

load resistance increases (Fig. 2. 9). Hie capacitors can't lose 
much of their charge when the load resistance is hi?^! .So the 
higher the load resistance , lower will be the peak to peak 
voltage ripple . 



3. 1.5. 4 Average output current follows the same pattern ais that 
of average output voltage . 

3. 1.5. 5 Strength of the source Current Inpulse 

This is directly proportional to the the peak to peak 
output voltage ripple and is given by C.AVo . Hence the effect of 
increased load resistance is same as that on peak to peak output 
voltage ripple . 

3. 1.5. 6 When the load terminals are open circuited tl^ following 
happens : Assuming we started with zero initial conditions ,at wt 
= 0* , there will be an impulse of current throu^ diode Dt and as 
a result Capacitor 1 will be charged to +E» immediately . Tire 
current through diode Di goes to zero after the capacitor 1 has 
been charged to peak value for there is no parallel path across 
capacitor 1 .Capacitor 1 S 2 remain charged with +Bb and 0 
respectively in the rest of the cycle from 0 to n . 

At (i)t = x*, Dz turns on an capacitor 2 is also charged to 
+Ei. owing to a negative source current iitpulse .The current 
throu^ Dz falls to zero immediately after the capacitor 2 has 
been charged . Capacitors 1 & 2 are thus dharged to 1.0 p.u. each 
and output remains at 2.0 p.u. for all values of time after ot = 
K^.^lre source current is also zero after ut = n*. 
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3.2 VOLTft® DOUBLING RECTIFIER FED El 

NOW IDERL SQURRE WRVE VOLTRGE SOURCE 

The doubling rectifier ( Fig. 3.7 ) fed by non ideal 
square wave voltage source has three nodes of cq?eration . Tlje non 
ideal square wave voltage source has resistance , R* and 
reactance, uLa . Tlie nodes are characterized by patterns of diode 
conduction . The nodes of operation in steady state are 

Mode 1 :Di ON , D 2 OFF, a< at i p ,a =0 

Mode 2 :Di OFF, Dz CffF, p < at ^ n and n + < at 2n 

Mode 3 :Di OTF, 02 0N,x<<otj:x+^ 

The description of each node is as follows : 

MCH% 1 : The equivalent circuit for this node is shown in Fig. 3. 9. 
In this node the diode Di starts conduction at at = 0 as the 
supply voltage is larger than capacitor 1 voltage . At ut = 0, a 
sinusoidal pulse of current flows through diode Di as the 
voltage across capacitor 1 is less than Em -Voltage across 
capacitor now jvnps to a value larger than Em (decided by R“L“C 
series coirbination) from a value less than Em as a result of this 
pulse of current. Diode Di turns off as soon as the source current 
becomes zero. The circuit enters Mode 2 at wt = p. The voltage 
across capacitor 2 falls- 



Fig. 3. 7 
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MODE 2 : The equivalent circuit for this mode is shown in 

Fig. 3. 10. The circuit enters mode 2 at wt = ^. In this mode both the 
diodes are off as they are reverse biased. Voltage across both the 
capacitors fall as they svq?ply energy to the load. 

MCOE 3 : The equivalent circuit for this mode is shown in 

Pig. 3. 11. The circuit enters mode 3 at ut = s as the 
negative of the sxjpply voltage is larger than the voltage across 
capacitor 2 or in other words diode Dz gets forw 2 u:xi biased. Hie 
voltage across capacitor 2 rises to a value larger than E* due 
to a negative source current pulse at ut = jt. Hie voltage across 
capacitor 1 falls in this mode. 

Hie circuit again enters mode 2 at ut = n+ p as the source 
current becomes zero. Hie diodes Di and D 2 remain off in this mode 
and the load is si^jplied by the stored energy in the capacitors. 

At ot = 2n, the mode 1 starts again and the same process repeats. 

Mode transition diagram in steady state is shown in Fig. 3. 8. 

3.2.1 DIFFERENTIAL BCPATIO® OF THE DODBLING RECTIFIER FOR MODES 
1, 2 AND 3 

In this section the differential equations describing the 
rectifier in different modes have been written . The derivatives 
of voltage across capacitors 1 S 2 and source current have been 
used in simulating the behaviour of the doubling rectifier . Hie 
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three variables of interest are Vci , Vc 2 and i, . The reference 
directicjn of capacitor voltage and current and source current are 
shovm in the Fig. 3.7 . The equations have been normalised as per 
the normalisation scheme given in the beginning of the chapter . 
Hie diodes have been assumed to be ideal . The pattern of modes in 
steady state is 1, 2, 3, 2, 1, 2 ... . The mode transition 

diagram is given in Fig. 3.8 . 

Mode 1 (Fig. 3.4) : In this mode Di is on and Da is off . Source 
current is positive . Hie circuit enters Mode 1 at wt = 0 . The 
three differential equations aure 

B, = i..R» + L*.d(i.)/dt + Vci [3.22] 

(Vci + Vc2) / R = -C.d(Vc2)/dt [3.23] 

source current, is = C.d(Vci)/dt + ( Vci + Vc2 )/R [3.24] 

Normalising the above equations , we have from equaticai 3.22 

1.0 = is.Rs/E« + (L«/Bs).[d(is)/dt] + Vci/Bs 

1.0 = (i8.Rs.<oC)/(B».wC) + (Ls.u^C)/(&.uC).[d(i«)/d(ut)] + Vci 

1.0 = isn.Rsn + (Ci)i«)n.^ d{ut)”] ^ [3.25] 

From equation 3.22 

( Vci + Vc2 )/El. = - (RC/&).d(Vc2)/dt 

Vein + Vc2n = ” (RC.«)/(B») •d(Vc2)/d(«t) 


.. B r dVc2n' 

Vein + Vc2n = - [ 


[3.26] 




Fig. 3. 9 Equivalent circuit in Mode 1 
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From equation 3.24 

i./(B«.wC) = C/(E;,.wC).d(Vci)/dt + (Vci + Vc2)/(R.B« .(oC) 

■ [-^Sr"] ♦ —- ’"fa [3-27] 

In all the equations derived above suffix n represents normalised 
variables e.g. («L*)n = o^.La.C r^resents normalised source 
reactance . 

Equations 3.22, 3.23, 3.24, 3.25, 3.26 and 3,27 describe the 
operation of voltage doubling rectifier in Mode 1 . 

MODE 2 (Fig. 3.10) : In this node both the diodes Di and D 2 are 
off . Hie circuit is characterised by source current , i» = 0 .The 
circuit enters mode 2 at ot = ^ vihafi source current becomes zero . 
The three differential equations are 


source current ,i» = 0 [3.28] 

(Vci + Vc2) / R = -C.d(Vc2)/dt [3.29] 

(Vcl + Vc2) / R = -C.d(Vci)/dt [3.30] 

Normalised equations can be written as follows 

ian ~ 0 [3.31] 


Equations 3.29 and 3.30 are of the kind of equation 3.23 , hence 
the nontalised equations can be directly written as 


_ r dVcln 

Vein + Vc2n = * 


[3.32] 
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Vein + Vc2n = - Rn. [3.33] 

Equaticjns 3.28 , 3.29 , 3.30 , 3.31, 3.32 and 3.33 describe the 
operation of the doiibling rectifier in node 2 . 

MODE 3 (Fig. 3.11) : In this node diode Di is off and diode Da is 
on . The circuit enters node 3 at ut = n when negative of the 
sv 5 )ply voltage becanes greater than voltage across capacitor 2. 
Hie source current is negative . The three differential equations 
of the circuit are derived below 

- EV, = i«.IU + U.d(i«)/dt - Vc2 [3.34] 

(Vci + Vc2) / R = - ad(Vci)/dt [3.35] 

source current, i* = - C.d(Vc 2 )/dt - ( Vci + Vc2 )/R [3.36] 

Normalisation of the equations can be done as earlier and the 
resulting equations in per unit are given below 


“ 1.0 — ian.Rsn + (oLs)n 


[3.37] 

B r dVcln] 

Vein + Vc2n = ' 


[3.38] 

. r d(Vcl)n'j j. 

- [ d((0€) J 

(V«t« + Vc2n) 

[3.39] 

Rn 


Equations 3.34, 3.35, 3.36, 3.37, 3.38 and 3.39 describe the 
cqperation of the doubling rectifier in node 3 . 
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3.2.2 Hie method of ccjnnputation. of numerical average output 
voltage and peak to peak output voltage ripple from the data 
points (instantaneous values of Vci, Vc 2 and ia over 2jt radians 
in steady state) is similar to that given in sec. 2.2.2. 

3.2.3 DISOJSSIOW OF SDfJLftTIOK RESOLTS In this section the 
general behaviour of the doubling rectifier has been studied . The 
transient bdiavi our, steady state bdtiaviour, average output voltage 
and peak to peak voltage ripple for different values of loaui 
resistances and source reactances have been calculated and shown 
in Figs. 3. 13, 3.14 and 3.15. These results are discussed below : 

3. 2. 3.1 L0iU3ED TORNSIEMT MID STERDT STATE BHfflVIOOR OF 1HE 
DODBLIMG RECTIFIER FED EY NON IDEAL SQOARE HAVE VOLTAGE SOORGE 

Hie circuit of Fig. 3.7 with a load resistance of 10.0 p.u. , 
source resistance of 0.01 p.u. and source inductive reactance of 
0.05 p.u. was simulated for studying its transient and steady 
state bdiaviour. 

Fig. 3.12 shows the simulated output voltage, voltage 
across capacitor 1, voltage across cpacitor 2, source current and 
iwp iif voltage. Hie circuit was excited from cold and behaviour 
recorded. 

At wt = 0*, the diode Dt turns on as the sv«)ply voltage is 
more than the voltage across capacitor 1. Capacitor 1 voltage at 
ut = 0* is zero as the circuit is excited from cold. Diode Da is 



off as it is reverse biased by the supply voltage. Differeace in 
supply voltage & the voltage across capacitor 1 at cot = 0*" caitseas 
a large sinusoidal current pulse. Biis pulse tends to becoaie sun 
impuse as the value of source reactance reduces . Due to this larg« 
pulse the voltage across capacitor 1 jumps to approadnatel y 1.^ 
p.u. 

During this period when diode Di is on diode D 2 is off , 
voltage across capacitor 2 increses in the negative directions, ats 
at the start of this interval the voltage across capacitor 2 vaas 
zero. 

The diode Di turns off soon as the source current becones 
zero and the circuit enters mode 2. In this mode of opecatiora, 
since the voltage across capacitor 1 is more than 1.0 p.u.^ EDi 
remains off and D 2 was already off. The two capacitors 
energy to the load and thus the output voltage falls. Hie voliage 
across capacitor 2 further increases in the negative direction arid 
touches -0.42 p.u. The voltage across capacitor 1 decreskses ircaii 
the level of 1.8 p.u. to approximately 1.5 p.u. 

Again at ut = ir , there appears a large negative saair«e 
current pulse. The amplitude of this pulse is much larger sas 
conpared to the pulse generated at wt = 0* , the reason kaernsr t3ie 
Icurge differencie in voltage across cpacitor 2 and tlM negative of 
the stpply. The voltage across capacitor 2 at wt = n actually 
to Em and causes this current pulse iduch raises the volt«ge 
across capacitor 2 to about 2.2 p.u. The rise in. the voltage is 
larger than the rise in the voltage across capacitor L at wt - 0. 
The source current soon crosses zero and diode D 2 turns of£. Di 




rol tage doubl ing 
roltage soxirce 



was off as it was reverse biased by the supply voltage. In this 
mode both the diodes are off and capacitors supply energy to the 
load. The circuit remains in this mode till wt = 2n 


MJLTIPLE CXMXJCTICti OF WE DKXES DURING WMSimT PERIOD 


Multiple conductiati of the diodes is a phenomenon in which 
the diode D 1 /D 2 conducts more than once in the positive/negative 
half cycles of the input sv^jply .Multiple conduction of the diode 
Di is noticed in the period from ot = 2it to 3ji. Multiple 
conduction of the diode D 2 is noticed in the period from ot = Sir 
to 4n. Multiple conduction occurs because of the fact that the 
diode D 1 /D 2 gets forward biased more tl^m once in the 
posit ive/negjative half cycles of the svqpply. 

In the period from wt = 2rc to Sir, the siqpply voltage is 
1.0 p.u. At the start of the interval there was a source current 
pulse and the voltage across cpacitor 1 rose to 1.3 p.u. Soon 
after the voltage rose to 1.3 p.u.. Diode Di went off and the 
voltage across capacitor 1 started falling and became equal to 1.0 
p.u. This forward biased the diode Di and it conducted again, 
second time, in the positive half cycle. 

Similar argiments apply for the explanation of the 
multiple conduction of Dz in the negative half circle from wt = 3k 
to 4rc. 

The circuit reaches steady state approximately after 
radians. In the staedy state, however, multiple conduction does 
not occur. In steady state the diodes conduction patterns is as 
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3. 2. 3. 4 Variation of the average output voltage vdth load 
resistance and source reactance 

Average output voltage vs. source reactance : 

The characteristics is shown in Fig. 3.14. From Fig. it is clear 
that the average output voltage falls both at very low values of 
wLa and at at relatively large values of wLa . The average output 
voltage is close to the ideal value of 2.0 p.u. for the range of 
uLa = 0.15 to 0.7. 

From the wave forms of the Fig. 3. 13 , the source current is 
seen to have sharp sinusoidal pulses at (aLs = 0.05 p.u. As uLa 
reduces towards ideal value of zero, these tend to become inpulses 
eis analysed in section in 3.1. Due to these sharp pulses of source 
current, the diodes Di and Dz turn off after conducting for very 
short intervals of time in the positive and negative half cycles 
respectively .Trtieref ore the capacitors si^jply energy to the load 
for most part of the cycle and thus output voltage falls for most 
part of the cycle. The average output voltage, hence, is low at 
low values of wLa . 

At high values of uLa ( > 0.7 p.u.) the peak source current 
reduces and due to large values of the source reactance the 
conduction period of the diodes in the positive and negative half 
cycles increases. However, the average output voltage falls as the 
inductive circuit is not able to transfer charge guickly. 

Average output voltage vs. load resistance .* 

The average output voltage rises as the load resistance is 
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increased. With the increased load resistance, the capacitors do 
not lose much of their stored energy and as a result the output 
voltage rises. 

3. 2. 3. 5 Variation of the Peak to Peak Output Voltage Ripple vd.th 
load resistance and source reactance 

Peak to Peak Output Voltage ripple vs. Source reactance: 

The characteristics is shown in Fig. 3.15. As e 3 q)lained in 
sec. 3. 2. 3. 4 above, the ripple voltage reduces as the source 
reactance is increased. At low values of the wLs, the diodes 
conduct large current pulses for shorter durations. These large 
current pulses cause sharp voltage rises across capacitors 1 & 2. 
For most part of the cycle diodes are off and capacitors discharge 
from their peak levels causing a large peak to peak output voltage 
ripple. 

At large values of source reactances, the conduction period 
of the diodes increase and the ripple in output voltage reduces. 

Peak to Peak output voltage ripple vs. load resistance: 

The peak to peak output voltage ripple reduces as the load 
resistance is increased. This is so because the capacitors 
voltage fall is less at higher values of the load resistances. 
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CJiPACITCR DIOCE MULTIPLIER EBD BY CURRENT SOURCES 

In this chapter, capacitor-diode voltage multiplier fed by square 
vrave and sine wave current sources has been stiKlied . linalysis is 
based on the differential equations modelling of the multiplier in 
its different modes of operation .Differential equations have been 
programmed and simulation results presented. 

4.1 CRPACITOR-DIOIE MULTIPLIER FED BY M 

IDEAL SINE WAVE AND SCARE WAVE CURRENT SOURCES 


The multiplier (Fig. 4. la and 4.1b) has two modes of operation 
characterised by the patterns of diodes conduction shown in Fig. 
4.2 . The modes are 

Mode 1 : Di on , D 2 off 0 < (ot i n 
Mode 2 : Da on , Di off 0 < «t s 2n 

HO PE 1 : In this mode ,the diode Di conduct sand Da remains off 
as the sr 5 >ply current is positive . The positive source current 
finds the path through capacitor 1 S 2 and the load resistance as 
shown in the Fig. 4.3 & 4.9 . The mode exists from ut = 0 to jt 
i.e. during the period when source current is positive . 

hO)E 2 : In this mode ,the diode Da conducts and Di remains off as 
the source current is negative. The negative source current finds 
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IDEAL SINE 
WAVE CURREViT 





Fig. 4. 1(a) Capacitor Diode Voltage Multiplier fed by an ideal 


sine wave current source. 



ideal 

Wave -SOORCE 


Fig. 4.1(b) Capacitor Diode Voltage Multiplier fed by aa ideal 
square wave current source. 


•+ 
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path through Dz and flows through capacitors 1 S 2 and load 
resistance as shown in Fig. 4.4 and 4.10 . The mode exists frcm o)t 
= n to 2r[ i.e. during the period when source current is negative . 
Cyclic changes in the modes of operation is shown in the Fig. 4. 2 . 
The patterns of modes is 1 , 2 , 1 and so on . 

4.2 DIFEBRSWTIIUL. BgOATICKS FOR THE CRERCITCR - DlCXa: MDLTIPLIHl 
FED BY IDEAL SINE WAVE CURRHTO SOURCE 


In this section the differential equations describing the 
capacitor diode multiplier in its modes of operation have been 
derived . The equations have been converted in p.u. form as per 
the nornalisation scheme given below : 

Normalisaticn Scheme : The applied sine wave current source is 
i«isin((i)t) . Anplitrjde im is taken as the base current . Ease 
inpedance (similar to the normalisation scheme chosen in earlier 
chapters) is 1/wC . Hence the base voltage is i«./uC . Normalised 
variables have been represented by suffix n . 


Mode 1 (Fig. 4.3) : This mode exists from 0 to n . Diode Di is on 
and D2 is off . The diffeential equations are 

i».sin(G)t) = C.d(Vci)/dt + (Vci + Vcz)/ R (4.1) 

- C,dj[Vc2)/dt = (Vci + Vc2)/ R ^^•2) 


Normalising the above equations as per the scheme given above 

sin(wt) = («C/i»).d(Vci)/d(ot) + {(Vci + Vc2).toC}/(R.i..oC) 


.(«t) = [.?!“] 
L d(«t)J 


+ (Vein + Vc2n)/Rn 


(4.3) 







READ INITIAL CONI^ITIOKJS / 
and load J^eSlST>^NiC6 / 

I ^ 

a)t = 0 INIT/ALI^E time CoUsllSK 


Canrpute and store the values 
of Vci , Vc2, is and Vo 
fron eq. 4. 3(4. 12), 4. 4(4. 13) 
for SINE (SQUME) wave current source 
with current, initial conditions 


Is 

.ut = re ?, 



+ tstep 


Compute and store the values 
of Vci, Vc2, is and Vo 
from eq. 4.7(4.16) ,4.8(4.17) 
for SINE (SQUARE) wave current source 
with current initial conditions 


ut + tstep 


Is 

ut = 2re ? 


/o(n+l)Tp - VoCnJTp = £ ? 

e = 0.01 / 


Using stored values of Vci, Vc2, is and Vo 
Print Average output voltage and j 

Peak to Peak voltage Ripple / 


Flow chart for Simulation of the capacitor diode 
multiplier fed by sine and square wave current source 


Fig. 4. 5 






- (c.u)/i«. r^fl = {(Vci + Vc2).(0C}/{R.i*.0)C) 
L d(wt)J 


- (Vein + Vc2n)/Rn (4.4) 


Equations 4.1 , 4.2 , 4.3 and 4.4 describe the (^>eration of 
capacitor - diode nultiplier fed by sine wave current source in 
mode 1 . 


Mode 2 (Fig. 4.4) : This mode exists from ut = n to 2n .In this 
mode D 2 is on and Di is off . The two differential equations are 


im.sinCut) + C.d[(Vc 2 )/dt + (Vci + Vc2)/ R = 0 (4.5) 

- C.d(Vct)/dt = (Vci + Vc2)/ R (4.6) 

Normalising the above equaticais as per the scheme given above 


sin((i)t) + (wC/im) .d(Vc2)/d((i)t) + {(Vci + Vc 2 ) .(>)C}/(R.i» .wC) = 0 

sin((i)t) + r + (Vein + Vc2n)/Rn = 0 (4.7) 

L'd(ut)J 

and 


- [ - (Vein + Vc2„)/Rn (4-8) 

L d(wt>l 

Equations 4.5 , 4.6 , 4.7 and 4.8 describe the operation of the 
multiplier in Mode 2 . 


4.3 DIFFERaJTIAL BQOATIONS FOR TSE CRPACITOR - DIODE MJLTIPLIER 
BED BY AN IDEAL SQUARE WAVE aJRRHff SOURCE 


In this section 


differential equations describing the 



multiplier in its inodes of operation have been derived . The 
equations have been converted in p.u. form as per the 
normalisation scheme in which peak current i» has been taken as 
the base current . Other beise quantities are same as that taken in 
last section . Normalised variables are denoted by suffix n . 


MOr® 1 (Fig. 4.9) : This mode exists from 0 to n . Diode Di is on 
and D2 is off .The differential equations can be written as 
follows 

i. = C.d(Vci)/dt + (Vci + Vc2)/ R (4.10) 

- C,d(Vc2)/dt = (Vet + Vc2)/ R (4.11) 

Normalising the above equations as per the schente given above 


1.0 = ((oC/i-).d(Vci)/d((ot) + {(Voi + Vc2).G)C}/(R.i-.(oC) 

1.0 = f + (Vein + Ve2n)/Rn (4.12) 

1 d(<i)t)J 

- (C.u)/ii».f = {(Vet + Ve 2 ).( 0 C}/(R.im.( 0 C) 

L“d(«t)J 

_ r dVc2n1 ^ + Vc 2 n)/Rn (4.13) 

Ld(«t)J 

Equations 4.10 , 4.11 , 4.12 and 4.13 describe the operation of 
capacitor - diode mrltiplier fed by square »ave current source in 

mode 1 . 

Mode 2 (Fig. 4.10) : This mode exists from et = a to 2x .In this 
„rx.e Da is cm and D, is off . «>e two differential equations are 
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- in, + C.d(Vc2)/dt + (Vci + Vc2)/ R = 0 


(4.14) 


- C.d(Vci)/dt = (Vci + Vc2)/ R (4.15) 


Normalising the above ecjuations as per the scheme given 


cuuruvc; 


and 


-1.0 + ((oC/im) .d(Vc2)/d({.)t) + {{Vci + Vc2).(oC}/(R.ii».{oC) = 0 

-1-0 + [_?!!"] + (Vein + Vc2n)/Rn = 0 (4.16) 

t d((0t)J 


dVcln 

d((i)t}- 


(Vein + Ve2n)/Rn 


(4.17) 


Equations 4.14 , 4.15 , 4.16 and 4.17 describe the operation of 
the multiplier in Mode 2. 


4.4 


METHOD OF SCaOTIOK 


The differential equations obtained for Mode 1 and Mode 2 have 
been solved numerically on corrputer using D02BAF SOBRODTINE from 
NAG library . The values of voltages across capacitor 1 S 2 and 
source current have been computed at a sufficiently large no. o§ 
points over a period from 0 to it in mode 1 and n to 2n for mode 2 . 
The computed values (herein after reffered to as data points) have 
been stored appropriately to be used later for computing average 
output voltage and peak to peak output voltage ripple . The 
flowchart for the program is shown in Fig. 4.5 . The transient 
behaviour has been recorded from zero initial conditions and 
steady state is reached after a few cycles . 

4.4.1 METHOD CF CX3MPOTATIC3N OP AVE3i2V3E C30TPDT VOLTAGE AND PHVK TO 
PEAK C30TfOT VOLTAGE RIPPLE 



This is ccrrputed in the manner as follows: 

The data points are available every ’tstep’ radians , 
tstep is the step size . We assume the output voltage to remain 
ccanstant over 'tstep' radians .The smaller the value of 'tstep' 
/higher will be the corrputaticsaal accuracy . Hence the average 
output voltage in per unit is 

Voavgn = V ._(Vclni tVc2ni) 
i.l " 

k = no. of data points in a period of 2n radians . 

In the conplete process of transient and steady state 
there can be pk data points where p is a large no. Csily last k 
datapoints are retained for these are the datapoints in which 
steady state has been reached . The datapoints during the 
transient state are only needed for plotting .The maximan and 
minirrajm values of the output voltage are pidced and the difference 
is outputed as peak to peak output voltage ripple . 

4.5 DISOJSSION OF SDfJIATION RESULTS 

In this section the transient and steady state brfiaviour 
of the rectifier has been recorded . Later in the section is 
stTjdied the variation of output voltage and peak to peak output 
voltage ripple with load resistance . 

4.5.1 LORDED 'nRMSIENT WD) STERDT STATE BEHAVICXJR OF IHE 
MULTIPLIER FED BY AN IDEAL SINE WAVE CURRENT SCORCE 

Fig. 4.6 shows the sMUlated output voltage , voltage 
across capacitor 1 , capacitor 2 and applied sinusoidal input 
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Fig. 4. 6(a) steady state behavioxir of the capacitor diode 
multiplier fed by sine wave current source 



current source for the rrultiplier having load resistance of 10.0 
p.u.The circuit of Fig. 4.1a is excited from cold and the 
behaviour is recorded . The circuit has two modes (as already 
explained in the beginning of the chapter)of operaticai 
Mode 1 . The capacitor 1 charges and capacitor 2 discharges 

(actually it charges in the negative direction ) during this 
mode. The voltage across capacitor 1 & 2 is almost sinusoidal . 

Mode 2 : Tlie reverse is true for mode 2 i.e. capacitor 2 charges 
and capacitor 1 discharges .The voltage across capacitor 1 & 2 are 
almost sinusoidal . 

Slowly the capacitors 1 & 2 build seme charge in them and 
steady state is reached after a certain no. of cycles .The voltage 
across capacitors 1 & 2 has a D.C. componait and a 100 Hz ripple . 
The D.C. compeanents add to give the average output voltage .The 
a.c. cotrponents have a phase difference but do not cancel each 
other exactly as in the square wave case . 

4.5.2 LQRE8ED TRftNSIEMT M® OTERDT STATE BfflAVIOCJR CF IHE 
MOLTIPDIER FED BY M lEERL SCPARE WAVE CORREST SCORCE 

Fig. 4.11 shows the simulated output voltage , voltage across 
capacitor 1 , capacitor 2 and applied square wave input current 
source for the multiplier having load resistance of 10.0 p.u.The 
circuit of Fig. 4.1a is excited from cold and the behaviour is 
recorded . The circuit has two modes (as already explained in the 
beginning of the chapter) of operation . 

Mode 1 : In the first half cycle the capacitor 1 charges as the 
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Fig. 4. 11 Loaded Traiisient and Steady State behavioiu: for the 
capacitor diode multiplier fed by square wave current 


source 




jupply current is positive . The capacitor 2 charges linearly in 
:he negative direction . At the end of n radians ,the circuit 
transits to mode 2 owing to the fact that source current changes 
Its direction . 

'lode 2 • The capacitor 2 voltage starts rising from the negative 
value from where it was left when the circuit transited to mode 2. 
It becanes positive soon . The voltage across capacitor 1 falls 
linearly from where it was left at wt = ji when the circuit 
transited to mode 2 . This process of charging and discharging of 
capacitors 1 St 2 respectively in mode 1 and the reverse in mode 2 
continues till the steady state is reached . In steady state the 
a.c. coni>onents of capacitors 1 S 2 voltages are 180° out of phase, 
llierefore the output voltage is nearly ripple free .The average 
ouput voltage is approximately 5.0 p.u. 


OPERATION AS A CURREm HALVER 

From the average output voltage vs. load resistance 
characteristics of the rectifier , it is clear that the circuit 
works as a current halver. At a load resistance of 10.0 p.u., the 
average output voltage for the rectifier fed by sine wave and 
square wave current sources is 4.23 p.u. and 4.98 p.u. 
respectively. Hence the load current is 0.423 and 0.498 times the 
source current respectively. The results show that the load 
current is , infact, lesser than half the ii^ut current for the 
rectifier fed by sine wave current source. This rmltiplication 
factor increases with the increase in load resistance and becomes 
0.5 at a load resistance of about 1000.0 p.u. However, for the 
souare wave current source the current halver oP®i^ution is bette 



and is seen at very low values of load resistances too. 

4.5.3 VARIATIONS IN AVERAGE OUTPUT VOLTAGE AND PEAK TO PEAK OUTPUT 
VOLTAGE RIPPLE WITH VARIATION IN LOAD RESISTANCE 

4. 5. 3.1 AVERAGE CXJTPUT VOLTAGE vs. LOAD RESISTANCE 

The average output voltage increases cts the load 
resistance is increased. It happens so because the source current 
is constant and any change in the output load resistance causes 
more current through capacitors. This enables them acquire more 
charge and as a result raise their voltages higher. In all 
this amounts to higher average output voltages. This is true for 
both the cases. Fig. 4.7 shows this behaviour. 

4. 5. 3. 2 REAK- TO PERK OUTPUT VOLTAGE RIPPLE vs. LORD RESISTANCE 
Fig. 4.8 shows this brfiaviour. 

(a) Square wave : The peak to peak output voltage ripple 

reduces as the load resistance is increased. The peak to peak 
output voltage ripple in square wave case is very low( zero for 
all practical pxirposes ). This is due to the fact the riKPles in 
vol tages across capacitor 1 & 2 are out of phase by 180 and hence 
cancel each other. 

(b) Sine wave : The ripple in the output voltage reduces as 

the load resistance is increased. This is due to the fact that at 
high values of load resistance, the variations in capacitors' 
voltages are saall since the current through than is high as 
explained in sec. 4. 5. 3.1 above and the charge loss is iiuninal 
because of the high load resistance. 
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CHAPTER 5 
SPICE SIMULATION 


The circuit of Fig. 2.1 has been analysed in chapter 2 & 
chapter 3 for operation as a voltage doubling rectifier. This 
circuit can be modified to obtain an output voltage which is 
adjustable in the range of 0 to 2 p.u. These modifications are 
discussed in this chapter. 

Need for controlling the output voltage arises in 
applications where regulation against variations in load and input 
supply is required. Some examples include HV power supplies for 
photomultiplier tubes used in scintillating counters, flying spot 
scanners etc. 

5.1 METHODS OF CONTROLLING THE OUTPUT VOLTAGE 

The voltage doubling rectifier of Fig. 2.1 can be modified in the 
following manners : 

(1) . The usual a.c. 50 hz input source can be replaced by an 
inverter operating from a fixed d.c. voltage source. This inverter 
can be modulated to obtain control of the output voltage. The 
control of the voltage rests in the Inverter while the doubling 
rectifier is unaltered. This cnfiguration is shown in Fig. 5.1. 

(2) An alternative is to modify the circuit of Fig. 2.1 by 
replacing the diodes by switches. It is then possible to modulate 
the input voltage to control the d.c. output voltage. This 
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Fig. 5.1 Voltage Multiplier fed by an Inverter operating on a fixed 
D.C source. 
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Fig. 5. 2 Voltage Doubling Rectifier fed by a Square wave Voltage 
source of controllable pulse width. 
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,he input voltage to control the d.c. output voltage. This 
configuration is shown in Fig.5.3. The square wave voltage is 
)ften conveniently generated by an inverter as shown in Fig.5.3 
<hile in principle sine wave input can also be used. The problem 
3 f reduced conduction due to charge on capacitors makes a square 
>iave input preferable. 

(3 The performance of capacitor - diode doubling rectifier fed by 
square wave and sine wave curent source has been given in chapter 
4. The circuit configurations are given in Fig. 4.1(a) & (b). The 
input source can be modified to a Pulse Width Modulated current 
wave to obtain controlled output voltage as shown in Fig. 5. 4 & 5.5. 

6.2 SIMULATION ON SPICE 

The circuits (1) to (4) above have been simulated on SPICE version 
3e2 available on main frame. SPICE is a computer aided circuit 
analysis program used for predicting the behaviour of electric and 
electronic circuits. Circuits are described to SPICE by use of an 
input file, which lists each circuit element (resistor, inductor, 
capacitance, voltage and/or current sources and semiconductor 
devices) and indicates how each is connected using node numbers. 
In addition there are lines in the input file which designate the 
frequency of the sources, temperature, the types of analyses to be 
lone and how the results are to be presented. The program also has 
nodel lines which describe the model of the semiconductor devices 

aeing used. 

The programs for the circuit of Fig. 5,1, 5.2, 5.3 and 
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isted in Appendix 3. The simulation of main subsystems in SPICE 
las been given below : 

,.2.1 SIMULATION OF MOSFET INVERTER AND MODULATION IN IT 

Enhancement mode NMOS is used as a switch in the full bridge 
inverter as shown in the Fig. 5.2. The MOSFET model chosen has the 
•ollowing details 

Level ( Model Index ) = 1 ( standard schichman hedges ) 

VTO ( zero bias thresshold voltage ) = 0 
Kp ( Transconductance Parameter ) = 1 
/TO and Kp define the d.c. characteristics of the device. The 
full bridge Inverter consisting of 4 NMOS devices has four voltage 
sources connected as shown In fig. 5.2, one for each NMOS device, 
to switch on and off each NMOS. The switching of each NMOS is done 
in such a way as to generate square wave voltage of controllable 
pulse width. For switching on/off of the device, the square wave 
voltage source of the type given in sec. 5.2.2 is used. 

5.2.2 SIMULATION OF THE PWM VOLTAGE SOURCE 

PWM voltage source has been simulated by connecting six 
voltage sources in series. Three voltage sources give positive 
pulses in the positive half cycle and three voltage sources give 
negative pulses in the negative half cycle. The frequency of the 
voltage generated by each source is 50 hz. The voltage source is 

defined in the following manner 

Vin = ( VI V2 Id Tr Tf PW Tp ) 

Vin » square wave voltage source of 
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Vi = voltage at t = Td‘ 

V 2 = voltage at t = Td* 

Td = delay time 

Tr = rise time 

Tf = fall time 

pw = pulse width (time during which voltage 
level is V 2 ) 

Tp= Time period of the voltage source 
The width in any of the voltage source can be varied to change the 
output voltage generated by the doubling rectifier. 

5.2.3 SIMULATION OF THE PWM CURRENT SOURCE 

PWM current source has been simulated by connecting six 
current sources in parallel. Three current sources give positive 
pulses in the positive half cycle and three current sources give 
negative pulses in the negative half cycle. The frequency of the 
current wave generated by each source is 50 hz. The current source 
is defined in the following manner 

lin = t II 12 Td Tr Tf PW Tp ) 
where lin = any periodic square wave current 

source of contollable width 
It = current at t = Td 

I2 = current at t = Td 

Td = delay time 

Tr = rise time 

Tf = fall time 
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level is I 2 ) 

Tp- Time period of the current source 
The width in any of the current source can be varied to change the 
output voltage generated by the capacitor diode multiplier. 

5.3 SIMULATION AND PERFORMANCE OF THE 

CIRCUITS DESCRIBED IN SECTION 5.1 
In this section the simulation of these controlled output voltage 
systems has been described. The following configurations have been 
studied. 

(1) The configuration shown in Fig. 5.2 has been studied. The 
voltage source inverter has a single pulse of controlled width 0. 

The pulse width can be modified in the range 0 to n to 
effect the change In the average output voltage and peak to peak 
output voltage ripple. The circuit of Fig. 5.2 with a load 
resistance of 500 ohms and capacitance of 100 pF was run. The 
performance of this circuit is given in the Fig. 5. 6(a) & (b). Some 
conclusions can be drawn from the simulation results: 

At large values of the load resistance , the effect of change 
in the pulse width 6 in the range from a to where a is very 
small, on the average output voltage is minimal. In other words, 
the output voltage is insensitive to the variations in the pulse 
width in this range. The output voltage, however, is sensitive to 
the variations in the pulse width in the range from 0 to a . The 
value of a is very low for large resistances. The smaller the 
resistance, higher will be the value of angle a. The value of « is 



actually decided by the circuit elements Rs and Capacitor 1 or 
2. If Rs.c is coparable to the pulse width of the input square 
wave, charging and discharging of the capacitors 1 & 2 can be 
controlled and hence the average output voltage can be controlled. 
With the reduced pulse width, source current peaks are higher and 
hence the output voltage ripple as explained in the sec. 3. 1.4. The 
sensitivity of average output voltage to the variations in pulse 
width increases with the reduction in load resistance. This 
behaviour is clear from Fig. 4.6(b) where variation in average 
output voltage is plotted for two values of the load resistance. 

(2) The performance of the doubling rectifier fed by a PWM voltage 
source with a load resistance of 1000 ohms and capacitance of 100 m 
has been studied and shown in Fig. 5. 8. The PWM scheme uses 3 
pulses /180 degrees and the frequency is fixed at 50 hz. The 
generation of PWM voltage source using SPICES has been discussed 
in sec. 5.2. 

(3) The configuration shown in Fig. 5. 3 with a load resistance of 
5000 ohms and capacitance of 25.5MF was simulated. The MOS 
inverter is running at 25 hz . The switches inside the multiplier 
were used to generate a controllable single pulse in each half 
cycle. The simulation of MOS inverter using SPICE3 is given in the 
section 5.2. Performance is given in Fig. 5. 7. 

(4) The configuration shown in Fig. 5. 5 with a load resistance of 
120 ohms and capacitance of 200 mF was simulated. The PWM current 
source has three pulses in 180 degrees and is operating at 50 hz. 
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Performance is given in Fig. 5. 9. The average output voltage can be 
better controlled In this scheme as the voltage build up across 
the capacitors is directly proportional to the charge supplied to 
them. By controlling the width of the pulses, we can control the 
voltage rise and fall across capacitors and hence the average 
output voltage. The simulation of PWM current source has been given 


in the sec. 5. 2. 
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Beised tipon the aiialysis and sinulaticn carried out in 
this thesis, the following conclusions can be drawn; 

(1) The full wave voltage doubling rectifier fed by an ideal sine 
wave voltage source has been analysed and simulated. Simulation 
results show that the circuit works as an acceptable voltage 
dotibling rectifier(average output voltage >1.9p.u.) for the values 
of load resistances greater than approximately 100.0 p.u. The 
performance degrades fast as the load resistance is reduced below 
100.0 p.u. Peak to peak output voltage rijple increases with the 
reduction in load resistance. 

(2) The full wave doubling rectifier fed by a square wave voltage 
source has been analysed and simulated. Simulation results show 
that its performance is better in terms of the output voltage 
regulation and peak to peak output voltage ripple as compared to 
the doubling rectifier in (1) above. However, the source curent 
peaks are much higher in this case. e. Diodes are found to conduct 
impulse of currents at wt = 0, n, 2n:. . 

3 (a) The voltage doubling rectifier fed by nonideal sine 

wave and square wave voltage sources has been analysed and 
simulated. The simulation results show that on a particualar 
value of wLs, the average output voltage is maxinun.The voltage 



falls as we increeuse the value of wLs above it or reduce the value 
of {i)Ls below it. This is found to be true for both nonideal sine 
and square wave voltage sources. 

(b) The effect of uLs on peak to peak output voltage 
ripple is different in the two systems. For the rectifier fed by 
nonideal sine wave voltage source, the peak to peak output voltage 
rippple increauses as the sovirce reactance is increased. For the 
rectifier fed by non ideal square wave voltage source, the peak to 
peak output voltage ripple decreases as the source reactance is 
increased. 

(c) Hie variation in peak to peak output voltage ripple 
for the doubling rectifier fed by non ideal sine wave voltage 
source is quite ccnplex. For low values of wLs (<0.04 p.u.), the 
increse in load resistance reduces the peak to peak output voltage 
ripple. For values of oLs larger than 0.04 p.u., the increase in 
load resistance till a particular value reduces the peak to peak 
output voltage ripple and a further increase in the value of load 
resistance beyond this value causes an increase in the output 
voltage ripple. 

Same thing is true for the variation in average output 
voltage with variations in load resistance .i.e. , average output 
voltage increases with increase in load resistance initially till 
a particular value of load resistance. If the load resistance is 
increased further, the output voltage falls.. 

For the rectifier fed hy a non ideal square wave voltage 
source, it is sirrpler to understand the bdiaviour of ripple 
voltage against variations in load resistance. The peak to peak 
output voltage ripple reduces with the increase in the load 



resistance. 


(4) The transient behaviour of the doubling rectifier fed 
the two voltage sources is also different. The rectifier fed by 
non ideal square wave voltage source shows the phenomenon of 
rmiltiple conduction of diodes Di and Dz in positive and negative 
half cycles respectively. The doubling rectifier fed by sine wave 
voltage source /however , does not show any such phenomenon. 

(5) The frequency of the ripple in the output voltage keeps on 
reducing as the load resistance is increased. For the Values of Rl 
lesser than 30.61 p.u. , the output voltage ripple frequency is 
twice as that of the sr 5 )ply frequency. This phaiomenon is not 
chserved in the rectifier fed by non ideal square wave voltage 
source. 

(6) The doubling rectifier fed by square and sine wave 
current sources has been analysed and simulated. The simulation 
results show that both the circuits work as a current halver. 

Voltage multiplier fed by sine wave current sources takes 
a large no. of cycles in reaching tte steady state . However , the 
multiplier fed by square wave voltage sotirce has a very small 
transient period. 

(7)Scme methods to control the output voltage have been 
suggested and simulated on SPICE. The doubling rectifier fed by a 
souare wave voltage source of controllable pulse width shoHs that 
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for large values of load resistance, the average output voltage is 
not much sensitive to the modulation of pulse width. The average 
output voltage does not change much in the range of the pulse 
width from a to jt where a is small and depends upon the 
product of R8,the equivalent source resistance, and the 

capacitance,C. However, in the range of the pulse width from 0 to 
a, the output voltage is very much sensitive to change in pulse 
width. 

As the load resistance is reduced the sensitivity of 
output voltage to the variations in pulse width increases. Output 
voltage ripple increcises as the pulse width is reduced in the 
pulse width range of 0 to 180 degrees. 

Modified voltage doubling rectifier and the rectifier fed 
by FWM voltage source also show the similar behaviour as 
explained above. The PWM current source gives much better control. 


SOGCTSTKWS PC» TOE FUTURE WORK 

(1) Mathematical analysis of the doubling rectifier' fed by a 
square wave voltage source of controllable pulse width can be done 
and expressions obtained for calculating the output voltage 
regulation and percentage ripple. 

(2) A better method of controlling the voltage douling rectifier 
has to be worked out for low values of load resistance. Detailed 
analysis of the effect of PWM voltage and current sources can be 


done. 
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APPENDIX NO. 1 

Derivations of the expressions for VC^(ot), and VC^<77+a) are 

given here. Equations 2.2, 2.5, 2.6, 2.11, 2.12, 2.16, 2.19 and 

2.20 are required for above derivations and hence are being 
reproduced here. 

In Mode 1 from o’jt = a. to 6jt = fiy we have 

VC, = E Sin oijt 2.2 

1 m 


VC. 


m 


2 2 2 
1+:VJ^R C‘ 


Sin iot + exp{-’t/RC) 


E viRC 
m 


2 2 2 
1+63 R C 


cos 03 1 


2.5 


K. 


' {'"=2 


E sin(oi) E 03 RCcos (ot )*. 

, , ffl m I 

ia) + — exp 


2 2 2 
1+03‘^R C*^ 


2^2:2 


1 +03 R C 


2.6 


Let 1 


4-^ and U 


E 03RC 
m 




1+03''^R'^C’^ 


then 


VC,. = 


sin 03t + K^exp C-t/RC) + cos 03t 


2 . 53 


K, = CVC,^(ot) + I, sin(ot) - i,, cosCaH 
X ^ 


2,63. 


In Mode 2 from 63t = ft to rr+a, we have 


VC., = E sin{/?) exp 
1 m 


[iis 


2.11 


^^2 = 1RC^3 


(il; 


2.12 


In mode 3 from ast = «+n to /?+?z, we have 


2, 16 


VC = - E sin 6jt 
Z m 


2^2 


VC^ ~ l^Esin iot ““ 60 R C cos 60t} + K^exp (~t/RC) 


= VC (rz+6?. ) 1 1 s in (6^ ) ~6 J^R'^C‘^cos (6?. ) ] exp — ;:^1 


In mode 2 again from 63t = 7Z+/? to 2n +ot, we have 


VCi = VC^<fZ+/?) exp ( 63 t- ( 71 '+/?) )J 


VC^ = Emsin(/?) exp (6it- ) )J 


From 2*51, we have 


VC^<27T+a) = Em sin(/‘?) exp 

= VC^(a) 


(ilr 


From 2.12, we have 


VC^ {«+/?) = VC^ (/?) exp <?J+a- 


(ilr- <"■*»-'’'] 


From 2.16, 


VC_(ot+n) = E sin{<a) 
2 m 


Hence 




E sin{a) exp 
m 


— (fT+a+/?M 

[rCoj - J 


From Eq. 2.11 


WC^in+a) = Em exp ^^^^- 


2. 19 

2. 20 

2. 50 

2.51 

2. 52 

2.53 

2.54 

2.55 

2, 56 


Equation 2,52, 2.55 and 2.56. give the expressions for 

VC (fi) and VC,{?T+ffl) respectively. 

.2 1 
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D02BAF - NAG FORTRAN Library Routine Document 

NOTI using this routine, please read the appropriate implementation document to check the interpretation of 

ttalictsea terms and other implementation-dependent details. The routine name may be precision -dependent. 


1, Purpose 

[X)2BAF integrates a system of first-order ordinary differential equations over a range with suitable 
initial conditions, using a Runge-Kutta-Merson method. 

2, Specification 

SUBROUTINE D02BAF (X. XEND, N. Y. TOL. FCN, W. IFAIL) 

C INTEGER N. IFAIL 

C realX, XEND, Y(N), TOL, WCN,7) 

C EXTERNAL FCN 


3. Description 

The routine integrates a system of ordinary 
differential equations 

Y/ = F.(T,YpY2....,Yj^), i = 1,2, ...N, 

from T = X to T = XEND using a Merson 
form of the Rungc-Kutta method. The system is 
defined by a subroutine FCN supplied by the 
user, which evaluates Fj in terms of T and 
Yj ,Y 2 ,...,Yj^ (see Section 5), and the values of 
Yj ,Y 2 ,....Yj^ must be given at T == X. 

The accuracy of the integration is controlled by 
the parameter TOL. 

For a description of Runge-Kutta methods and 
their practical implementation see [1]. 

4 . References 

(1] HALL, G. and WATT, J.M. (eds). 

Modern Numerical Methods for Ordinary 
Differential Equations, p. 59. 

Clarendon Press, Oxford, 1976. 

5. Parameters 
X - real. 

Before entry, X must be set to the initial value 
of the independent variable T. 

On exit, it contains XEND, unless an error has 
occurred, when it contains the value of T at the 
error. 

XEND - real. 

On entry, XEND must specify the final value 
of the independent variable. If XEND < X on 
entry, integration will proceed in the negative 
direction. 

Unchanged on exit. 


N - INTEGER. 

On entry, N must specify the number of 
differential equations. 

Unchanged on exit. 

Y - real array of DIMENSION at least (N). 

Before entry, Y(1),Y(2),...,Y(N) must contain 
the initial values of the solution Yj ,Y 2 v-»Yj^ . 

On exit, Y(1),Y(2),...,Y(N) contain the 
computed values of the solution at the final 
value of T, 

TOL - real. 

Before entry, TOL must be set to a positive 
tolerance for controlling the error in the 
integration. 

The routine D02BAF has been designed so that 
for most problems a reduction in TOL leads to 
an approximately proportional reduction in the 
error in the solution at XEND, However, the 
actual relation between TOL and the accuracy 
achieved cannot be guaranteed. The user is 
strongly recommended to call D02BAF with 
more than one value for TOL and to compare 
the results obtained to estimate their accuracy. 
In the absence of any prior knowledge, the user 
might compare the results obtained by calling 
D02BAF with TOL = 10.0“^ and 
TOL = 10.0 if P correct decimal digits 
in the solution are required. 

TOL is normally unchanged on exit. However, 
if the range X to XEND is so short that a small 
change in TOL is unlikely to make any change 
in the computed solution then, on return, TOL 
has its sign changed. This should be treated as 
a warning that the computed solution is likely 
to be more accurate than would be produced by 
using the same value of TOL on a longer range. 
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f C'N SUBROUTINE, supplied by the user. 
FCN must evaluate the functions (i.c. the 
derivatives Yj' ) for given values of its 

arguments T,yj . 

Its specification is: 

SUBROUTINE FCN(T. Y.F) 
rea/ T.Y(n),F(n) 

where n is the actual value of N in the call of 
D02BAF. 

T - real. 

On entry, T specifies the value of the 
argument T. 

Its value must not be changed. 

Y real array of DIMENSION (n). 

On entry, Y(I) contains the value of the 
argument Yj , for I = l,2,...,n. 

These values must not be changed. 

F -> real array of DIMENSION (n). 

On exit, F(I) must contain the value of Fj, 

I = 

FCN must be declared as EXTERNAL in the 
(sub)program from which D02BAF is called. 

W - real array of DIMENSION (N,p), where 
p> 7. 

Used as working space. 

IFAIL ~ INTEGER. 

On entry, IFAIL must be set to 0 or 1. For 
users not familiar with this parameter 
(described in Chapter POl) the recommended 
value is 0. 

Unless the routine detects an error (see next 
section), IFAIL contains 0 on exit. 

6. Error Indicators and Warnings 

Errors detected by the routine:- 

IFAIL = 1 

On entry, TOL :< 0.0 
or N :< 0. 

(The latter error will cause a program 
Weakdown with some compilers.) 

IFAIL = 2 

With the given value of TOL, no further 
progress can be made across the integration 
range from the current point T = X, or the 
dependence of the error on TOL would be 


lost if further progress across the integration 
range were attempted (see Section 11 for a 
di.scussion of this error exit). The components 
Y(1),Y(2),...,Y(N) contain the computed 
values of the solution for the current point 
T = X. 

IFAIL - 3 

TOL is too small for the routine to start the 
integration (see Section II), X and 
Y(I),Y(2),...,Y(N) retain their initial values. 

IFAIL = 4 

A serious error has occurred in an internal 
call to D02PAF. Check all subroutine calls 
and array dimensions. Seek expert help. 

7. Auxiliary Routines 

Details are distributed to sites in 
machine-readable form. 

8. Timing 

This depends on the complexity and 
mathematical properties of the system of 
differential equations defined by FCN, on the 
length of the range and on the tolerance. There is 
also a small overhead of the form A + B X N, 
where A and B are machine-dependent 
computing times. 

9. Storage 

The storage required by internally declared 
arrays including those of auxiliary routines, is 28 
real elements. 

10. Accuracy 

The accuracy depends on TOL, on the 
mathematical properties of the differential 
system, on the length of the range of integration 
and on the method. It can be controlled by 
varying TOL but the approximate 
proportionality of the error to TOL holds only 
for a restricted range of values of TOL. For TOL 
too large, the underlying theory may break down 
and the result of varying TOL may be 
unpredictable. For TOL too small, rounding 
errors may affect the solution significantly and 
an error exit with IFAIL = 2 or IFAIL = 3 is 
possible. 

The user who requires a more reliable estimate of 
the accuracy achieved than can be obtained by 
varying TOL, is recommended to call the routine 
D02BDF where both the solution and a global 
error estimate are computed. 
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11. Further Comments 

If Ihc routine fails with IFAIL = 3, then it 
could be called again with a larger value of TOL 
(if this has not already been tried). If the 
accuracy requested is needed and cannot be 
obtained with this routine then the system may 
be very stiff (see below) or so badly scaled that it 
cannot be solved to the required accuracy. 

If the routine fails with IFAIL = 2, it is 
probable that it has been called with a value of 
TOL which is so small that a solution cannot be 
obtained on the range X to XEND. This can 
happen for well-behaved systems and very small 
values of TOL. The user should, however, 
consider whether there is a more fundamental 
difficulty. For example, 

(i) in the region of a singularity (infinite 
value) of the solution, the routine will 
usually stop with IFAIL = 2, unless 
overflow occurs first. If overflow occurs 
using D02BAF, routine D02PAF can be 
used instead to trap the increasing 
solution before overflow occurs. In any 
case, numerical integration cannot be 
continued through a singularity, and 
analytical treatment should be 
considered. 

(ii) for ‘stiff equations, where the solution 
contaiins rapidly decaying components, 
the routine will use very small steps in T 
(internally to D02BAF) to preserve 


stability. This will exhibit itself by 
making the computing time excessively 
long, or occasionally by an exit with 
IFAIL = 2. Merson’s method is not 
efficient in such cases, and the user 
should try the Gear method D02EAF. 

To determine whether a problem is stiff, the 
routine D02BDF may be used. 

Users with problems for which D02BAF is not 
sufficiently general should consider the routines 
D02BBF, D02BHF or D02PAF. Routine 
D02BBF can be used when output is required at 
points outside the range X to XEND (for 
example, for graph-plotting purposes) or more 
general error control is requir^. Use of IX)2BBF 
should be computationally more efficient than 
repeated calls to D02BAF to achieve the same 
result. Routine D02BHF can be used to calculate 
where a function of the components 
Yj and their derivatives take a 

specified value. 

I>02PAF is a more general Merson routine with 
many facilities including more general error 
control options and several criteria for 
interrupting the calculations. 

12. Keywords 

Initial Value Problems, 

Ordinary Differential Equations, 
Runge-Kutta-Merson Method. 



PROGRAM NO. 1 


voltage doubler fed from a mos inverter 
vin 1 0 dc 50 

vinl 2 3 pulse(-2 10 Om O.lm 0.1m 10m 20m) 

vin2 4 5 pulse(-2 10 10m O.lm O.lm 10m 20m) 

vin3 7 0 pulse(-2 10 Om O.lm O.lm 10m 20m) 

vin4 6 0 pulse(-2 10 10m O.lm O.lm 10m 20m) 

mosl 1233 simple 

mos2 1455 simple 

mos3 3600 simple 

nK)s4 5 7 0 0 simple 

rs 3 9 1 

dl 9 10 d 

d2 11 9 d 

cl 10 5 lOOu 

c2 5 11 lOOu 

rl 10 11 1000 

*re 5 11 lOOmeg 

*rel 11 0 200meg 

*re2 4 0 700meg 

*re3 3 0 200meg 

*relO 10 0 200meg 

.opticjns nopage 

.model simple nmos(kp=l,vto=l,level=l) 

.model d d 

.tran Im 60m 

.plot tran v(10,ll) 

.end 


PROGRAM NO. 2 

voltage doubling rectifier fed by mos inverter 
*diodes in the doubling rectifier replaced by switches 
vin 10 0 dc 50 

vinl 20 9 pulse(“2 10 0 0.02m 0.02m 2Sn 50m) 
vin2 14 0 pulse(-2 10 25m 0.02m 0.02m 2&n 50m) 



vin4 15 0 pulse(-2 10 0 0.02in 0.02m 25m 50m) 

vgl 2 3 pulse(-2 10 5m 0.02m 0.02m 20m 50m) 

vg2 6 1 pulse(-2 10 30m 0.02m 0.02m 20m 50m) 

*mosfet switches 

mosl 10 12 9 9 simple 

mos2 9 14 0 0 simple 

mos3 10 13 7 7 simple 

inos4 7 15 0 0 simple 

m2 5 6 1 1 simple 

ml 1 2 3 3 simple 

cl 4 7 25. 5u 

♦reverse recovery diodes 

da 19 13 dl 

db 20 12 dl 

dp 9 10 dl 

dq 7 10 dl 

dr 0 9 dl 

ds 0 7 dl 

rsl 9 1 5 

rl 4 11 5000 

c2 7 11 25. 5u 

*mega ohm resistors provided 

*so that each node 

♦gets a DC path to ground 

rel 7 11 lOOmeg 

re2 9 0 lOOroeg 

re3 7 0 lOOmeg 

re4 1 0 lOOmeg 

re5 3 0 lOOmeg 

re6 5 0 lOOmeg 

re7 12 0 lOOmeg 

reS 19 0 900meg 

dl 11 5 dl 

d2 3 4 dl 

.tran Im 100m uic 
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.model sinple nmos(vto=l,kp=l,level=l) 
.model dl d 
.options nopage 
.plot tran v(4,ll) 

.plot tran v(9,7) 

.end 


PROGRAM NO. 3 

voltage doubler FVW MALYSIS 

■* definition of a PWM voltage source 

vinl 1 2 pulse(0 50 l.llm O.Olm 0.01m 1.11m 20m) 

vin2 2 3 pulse(0 50 3.3au O.Olm O.Olm 3.33m 20m) 

vin3 3 4 pulse(0 50 7.77m O.Olm O.Olm l.llm 20m) 

vin4 4 5 pulse(0 -50 11.11m O.Olm O.Olm l.llm 20m) 

vin5 5 6 pulse(0 -50 13.33m O.Olm O.Olm 3.33m 20m) 

vin6 6 0 pulse(0 -50 17.77m O.Olm O.Olm l.llm 20m) 

♦source resistance 

rs 1 8 1 

♦diodes and capacitances 

dl 8 9 dl 

cl 9 0 lOOu 

c2 0 10 lOOu 

d2 10 8 dl 

♦load resistance 

rl 9 10 1000 

♦default model of the diode 
.model dl d 
.options nopage 
♦transient analysis 
.tran O.lms 80ms 

♦plot output voltage and input voltage 
.plot tran v(9,10) 

.plot tran v(l) 

.end 
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PROGRAM NO. 4 

voltage doubler EVM CURREOT SOURCE ANALYSIS 
*doubling rectifier fed by PHM current vrave 

* definition of the current source 

iinl 1 0 pulse(0 1 l.llm O.Olm 0.01m l.llm 20m) 

iin2 1 0 pulse(0 1 3.33m O.Olm O.Olm 3.33m 20m) 

iin3 1 0 pulse(0 1 7.77m O.Olm O.Olm l.llm 20m) 

iin4 1 0 pulse(0 -1 11.11m O.Olm O.Olm l.llm 20m) 

iinS 1 0 pulse(0 -1 13. 3M O.Olm O.Olm 3.3!^ 20m) 

iin6 1 0 pulse(0 -1 17.77m O.Olm O.Olm l.llm 20m) 

♦source resistance 
rs 1 7 1 

* zero volts voltage source 
vdo 780 

♦diodes and capacitors 
dl 8 9 dl 
cl 9 0 200u 
c2 0 10 200u 
d2 10 8 dl 

♦default model of the diode 

rl 9 10 120 

.model dl d 

.options nopage 

.tran O.Sm 140m 

.plot tran v(9,10) 

.plot tran i(vdo) 

.end 
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